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Extensions of the MSSM often predict the existence of new fermions and their scalar superpartners 
which are vectorlike with respect to the standard model gauge group but may be chiral under 
additional gauge factors. In this paper we explore the production and decay of an important example, 
i.e., a heavy isosinglet charge —1/3 quark and its scalar partner, using the charge assignments of a 
27-plet of Es for illustration. We emphasize that, depending on the symmetries of the low energy 
theory, such exotic particles may decay by the mixing of the fermion with the d, s, or b quarks; may 
decay by leptoquark or diquark couplings (which may nevertheless preserve a form of ii-parity); or 
may be stable with respect to renormalizable couplings but decay by higher-dimension operators on 
cosmological times scales. We discuss the latter two possibilities in detail for various assumptions 
concerning the relative masses of the exotic fermions, scalars, and the lightest neutralino, and 
emphasize the necessity of considering the collider signatures in conjunction with the normal MSSM 
processes. Existing and projected constraints from colliders, indirect experiments, proton decay, 
and big bang nucleosynthesis are considered. 

PACS numbers: 12.60. Jv,04.65.+e,14.80.Ly,95.35.+d 
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I. INTRODUCTION 

As the beginning of data-taking at the CERN Large 
Hadron Collider (LHC) rapidly approaches, there has 
been an increase in research into both the possible new 
physics opportunities as well as the challenges facing ex- 
perimenter and theorist alike. The focus of attention has 
thus far been squarely on new physics of a supersymmct- 
ric nature, in particular on the minimal supersymmetric 
extension of the Standard Model (MSSM). This is un- 
doubtedly a well-justified approach, but in our prepara- 
tory studies of opportunities and challenges at the LHC 
it is wise to sometimes take a broader view of what the 
TeV scale may reveal. There has recently been a grow- 
ing appreciation 

0, H H II IS 1 of the possible difficul- 
ties in connecting data to underlying theories, focusing 
on extracting parameters of the MSSM Lagrangian or 
distinguishing the MSSM from some other paradigm of 
new particle physics. In this paper we wish to consider 
another possibility. We will consider here the case in 
which the "usual" new physics signal from the states of 
the MSSM is combined with that from additional super- 
symmetric states accessible at the TeV scale. In particu- 
lar, we will consider the presence of additional vector- like 
pairs of SU (3) triplets which are singlets under the SU (2) 
of the Standard Model but carry non-vanishing hyper- 
charge (as well as possible additional U(l)' charges). 

Such states, which we will denote as D and D c in 
this work, are among the most well-motivated extensions 
of the MSSM particle content. For example the super- 
symmetric SU(5), SO(10) and Eq grand unified theories 



(GUTs), as well as non-supersymmetric Eq, predict the 
existence of such states 0, Q, though only in the Eq 
case (or in nonminimal versions of SU(5) and 5O(10)) 
can some of them have electroweak-scale masses without 
severely compromising the lifetime of the proton. Simple 
enlargements of the electroweak sector of the MSSM in- 
corporate these states, to saturate anomaly cancelation 
constraints in U(\)' models and/or to contribute to the 
radiative breaking of the extended symmetries 0, [ljl HH, 
El, [H, 13- extensions 1 may arise, for example, 

as variations on the NMSSM [l5[ for generating dynam- 
ical \x terms. Structures of these sorts arise routinely in 
top-down string constructio ns, part icularly those of the 
heterotic string 0, 03, [H EHM HH [H, HI l| - so 
much so, in fact, that the -EVversion of D and D c are of- 
ten referred to simply as "superstring-inspired" exotics. 2 
Yet more significant (for our purposes) than these moti- 
vations is the simple fact that the existence of additional 
vector-like pairs of SU(3) triplets is a very real logical 
possibility that will have a significant impact on the ex- 
perimental environment at a hadron collider such as the 
LHC. 3 



Such extensions are consistent with simple gauge coupling uni- 
fication if the D and D c occur in combination with additional 
vector-like Higgs or lepton doublet pairs, i.e., so that the new 
states have the quantum numbers of a 5 + 5* of 5(7(5). 
In explicit string constructions one finds that additional 
triplcts/antitriplets of SU(3) often have U(l)' charges that dis- 
tinguish these states from their GUT counterparts. 
In this regard the present work serves as a complement to similar 
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Consider, for example, the various logical possibilities 
for interaction between exotic D and D c states and the 
fields of the Standard Model. The simplest scenario, and 
the one that is usually studied, especially in the non- 
supersymmetric case, is one in which the exotic D and the 
Standard Model d are allowed (by the conserved quantum 
numbers) to mix, as are the conjugate fields D c and d c . 
Another scenario might forbid this sort of mixing (e.g., 
because of R-parity conservation) , but allow renormaliz- 
able interactions between the new D, D c and Standard 
Model fields in the superpotential. Depending on the 
quantum numbers of the exotic states a range of pos- 
sible production and decay mechanisms would then be 
operative, including interactions that are leptoquark or 
diquark in nature. Still another possibility is that some 
new symmetry - or combination of symmetries - for- 
bids these renormalizable interactions but allows higher- 
order interactions in the superpotential. In this case the 
new exotics would be "quasi-stable" ; they would have a 
lifetime that implies stability on collider timescales (ap- 
proximately 100ns or longer) but might decay sufficiently 
quickly to avoid cosmological limits on such objects. Fi- 
nally, we mention the possibility that the exotic particle 
is absolutely stable or stable on the time scale of the age 
of the universe. This would imply cosmological difficul- 
ties and is also severely constrained by searches for heavy 
stable particles, e.g., in sea water (for a review, see [36|). 
However, these difficulties could possibly be avoided if 
the reheating temperature of the universe after inflation 
was very low [37[ , especially if the mass is larger than a 
TeV [38|. 

In this work we will therefore cover a variety of such 
phenomena: new signatures and new challenges for con- 
necting data to theories, cosmological issues of long-lived 
heavy 5/7(3)-charged states, and indirect constraints 
coming from limits on rare processes. The goal is to 
be as comprehensive as possible in enumerating all the 
logically distinct ways in which these exotic quarks and 
squarks can manifest themselves in Nature. In each case 
we wish to ask what are the constraints on the existence 
of these new states, and by what observational methods 
will we infer their existence? We would also like to know 
how the answers to these questions will depend on the 
free parameters of theories which contain such exotics. 
To this end we will need to consider an array of possi- 
bilities instead of any one paradigm. Nevertheless, to be 
concrete, one needs to discuss phenomena in terms of a 
model or class of models. We will therefore choose rep- 
resentative ones when necessary, but we strive to treat 
everything as model-independently as is possible. We 
note that several specific models that fall into this class 
have already been studied at length, but other aspects 



which we cover have hardly been mentioned in the liter- 
ature. More complete references to existing studies will 
be provided in the appropriate sections. 

The rest of the paper is organized as follows. In Sec- 
tion HQ we will provide some background on iVbased 
models and explain the various ways that exotic isosin- 
glet quarks and squarks can arise. Much of this material 
will be a review of previous work. In Section [II CI we will 
provide five benchmark scenarios defined by the masses 
of the exotic fermion and its scalar superpartners. These 
benchmark cases are designed to cover a range of likely 
soft supersymmetry breaking mass scales and subsequent 
phenomenology. In Section III Dl we take a first look at 
production processes for these exotics at hadron collid- 
ers, and discuss the modifications we made to the PYTHIA 
event generator to handle these new states. Section IIIII 
summarizes the bounds on exotic SU (3)-charged isosin- 
glets arising from direct searches, rare flavor-changing 
processes and cosmology. The latter are constraining 
only for cases in which the exotic states are quasi-stable. 
In Section IIVI we look at the collider phenomenology of 
this quasi-stable case in great detail, and give the discov- 
ery reach for such states at the LHC. Section|V]is devoted 
to the collider phenomenology of scenarios in which the 
exotic decays promptly in the detector through renor- 
malizable interactions. For technical reasons we will here 
focus on the more tractable case in which the exotic has 
leptoquark-type couplings, reserving the diquark case for 
a separate work. Additional material is contained in two 
appendices. 



II. OUTLINE OF CASES 

In this section we will present some basic model con- 
cepts that will allow us to treat the phenomenology of 
exotic D and D c quarks in a semi-unified manner. The 
smallest extension of the MSSM capable of incorporat- 
ing all of the phenomena mentioned in Section |T] is mo- 
tivated by Eg GUT symmetries. Systems based on the 
Eq gauge group (or its various subgroups), with mat- 
ter states arising from the fundamental 27 representa- 
tion, have been in the past - and continue to be now - 
an appealing framework for organizing models of beyond 
the MSSM physics [U, M, El, E3, El 0, SI, [H, Sl- 
From within this framework many interesting limits and 
sub-models can be studied: quasi-stable exotics, exotics 
which mix with SM fermions, promptly decaying lepto- 
quarks, promptly decaying diquarks, and so forth. Let 
us emphasize that we are not interested in any particu- 
lar Eq model, nor do we even insist that all elements of 
grand unification are present. 4 We will also not con- 



treatments of new SU (2) doublets and new Higgs singlets [la . 

SllH,m,[I3,m[2l[il,li or other types of new P h y sics 

at the LHC. 



4 For example, if the D and D c states related to the Higgs doublets 
responsible for fermion mass generation are at the TeV scale (as 
is required if there is a TeV scale U(l)' gauge symmetry), then 



cern ourselves much with the additional richness that 
such models present (new Z'-bosons, new Higgs dou- 
blets, additional neutralinos and Higgs singlets, chal- 
lenges for neutrino mass generation and gauge coupling 
unification, etc.). We merely introduce some basic facts 
of _Eg-inspired models for the sake of coherence of pre- 
sentation and as an example of consistent, anomaly- free 
exotic particle quantum numbers that can be consistent 
with gauge unification. Many constructions which do not 
fit into the E% framework can give rise to the physics we 
will describe in Sections IIVI and [V] and our discussion of 
the consequences of exotic D and D c pairs will be quite 
general. 
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A. General Framework 

The Eq gauge group decomposes to the gauge group 
of the Standard Model as 

E 6 -> SO(10) x 17(1)^ 

-> SU(5) x U(l) x x U{1)^ (2.1) 
-> SU(3) x SU(2) x U(1) Y x U(l) x x U(l)^, 

where the designation of the particular U(l) combina- 
tions and U(l) x are traditional and chosen for 
convenience of notation. The above symmetry breaking 
can occur by a variety of physical means. For our pur- 
poses the exact mechanism is unimportant, though we 
might imagine a string-theoretic origin for our exotics, 
in which case a Wilson line breaking analogous to the 
Hosotani mechanism may be envisioned (48|. There are 
two additional U(l) factors beyond those of the Standard 
Model. One or two linear combinations of these factors 
may remain intact to very low energies, though both must 
be broken at the electroweak scale. We will assume the 
associated Z'-bosons are heavy enough to be irrelevant 
to the LHC phenomenology we wish to explore. 

Each fundamental representation 27 of E§ contains the 
Standard Model representations given in Table Q] The 
particle content of the 16 representation of 50(10) is 
augmented by a pair of Higgs doublets, a pair of ex- 
otic quarks D and D c and a singlet field S. This par- 
ticle content is anomaly-free by construction. Achiev- 
ing three generations of Standard Model fields there- 
fore implies three generations of Higgs fields, exotic 
triplcts/antitriplets and singlets. With the field content 
of Table H] it is well known that gauge coupling unifi- 
cation cannot be achieved consistent with the measured 
low-energy values of , a EM and sin 2 9w ■ This can be 
easily remedied by the introduction of additional fields 
in an anomaly- free manner [l2l [4^ . [Hoj ] . As this is im- 



the E(j relations between the Higgs and D, D c Yukawa couplings 
must not be respected in order to avoid rapid proton decay. For- 
tunately, string constructions often do not honor such relations. 



TABLE I: Decomposition of the fundamental of Ee un- 
der the Standard Model gauge group. The fields of 
the 27 i representation of Eq are here given in terms of their 
SU(3) x SU(2) representation, as well as their charges under 
the four (7(1) factors. (Qn is relevant to the quasi-stable sce- 
nario.) The index i = 1, 2, 3 represents a generation index. 
Some of the H Ut d pairs could be interpreted as exotic lepton 
doublets if they do not acquire expectation values. 

material to our purposes we will not consider the issue 
further. 

The purpose of introducing the field content of the 27 
representation was simply to motivate the form that su- 
perpotential interactions of D and D c might be allowed 
to take. The renormalizable superpotential for Eq is sim- 
ply given by W = Ayfc27j27 J -27/ c . However, we do not 
insist on full Eq invariance of the superpotential, but 
rather use this form to motivate the classes of allowed 
couplings. When decomposed into the appropriate com- 
ponents under (|2.1[) W then becomes 

W = x],(lu\n„ ■ A;,o,-/;//,, • x;,i.y,n, 

• A;; 1.,,/JI,, + \*SH«H d + Xl-SDiD* 
+W LQ + W UQ (2.2) 

where 

W LQ = \% k D lU c 3 ei + ^DtQjL* + \% k D % d c 3 v L k (2.3) 
^dq = X^DiQjQk + X^D^dt, (2.4) 

and the convention for numbering the interactions is 
taken from the review of Hewett and Rizzo [8| . We have 
restricted our attention here to one relevant generation of 
each. This could be achieved by appropriate assumptions 
concerning the Yukawa matrices and/or the vacuum ex- 
pectation values (vevs) of the singlet fields Si. We will 
usually restrict our attention to a single generation of D 
and D c quarks as well, but for now we retain the gener- 
ation label. 

If we were to demand invariance only under the Stan- 
dard Model gauge group then additional bilinear and 
trilincar terms (e.g., terms involving just the standard 
model fields which violate R-parity) could be added 
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to (|2.2p [5ll]. 5 If we require invariance under only the 
Standard Model plus one additional £7(1) factor then 
some subset of these additional terms may be allowed. 
If one additional U(l) factor arising from the original E§ 
is retained to low energies it is traditionally parameter- 
ized as 

Q' = Q x cos#e + sm9 E , (2.5) 

where the charges Q x and Q^p are those given in Ta- 
ble |U Any choice of 9e in (|2.5p allows all the terms 
in (|2.2p - (|2.4p . by construction. Higher-dimensional, non- 
renormalizable operators are also possible in the super- 
potential. Their presence or absence depends on which 
linear combination in (|2.5p . if any, is assumed to be 
present at low-energies. For the sake of concreteness, 
when necessary we will choose U(l)' charge assignments 
for these fields according to the [7(1)^ combination with 
6e = 2-7T — tan -1 \/5/3, or to the U(1)n combination 
with 9e = tan -1 yTE. 

B. Charge Assignments 

If both (|2.3p and (|2.4p are present simultaneously then 
it is impossible to assign an unambiguous B and L quan- 
tum number to D and D c - thus B and L are broken. 
In this case the exotic SU (3)-charged states will mediate 
rapid proton decay. We will therefore insist on sepa- 
rately conserved quantum numbers B and L and choose 
superpotential terms to allow definite B(D) and L(D) as- 
signments. This will always imply a trivially conserved 
R-parity quantum number using the standard definition 

i? p = (ll)3(S-i)+ 2s 

When only (|2.3p is present then one can assign the 
quantum numbers B(D) = 1/3 and L(D) = 1 so that 
R P {D) = — 1 and R P (D) = +1 and we can identify D 
as a leptoquark. With only (|2.4[) we have B(D) = —2/3 
and L(D) = and the same R p assignment; the state D 
is then a diquark. Note that in these two cases the D 
and D c are like H u and the scalar is the "standard" 
particle and the fermion is the "partner" . So the R-parity 
distinguishes 5's associated with the 16 of SO(10) from 
those coming from the 10's of SO(10). In this case the 
only renormalizable operators allowed are those of (|2.2p 
with (|2.3p or (|2.4p . All dimension- five operators involving 
the exotic D and D c also vanish in this case. 

If we instead insist on baryon and lepton number con- 
servation with B(D) = 1/3 and L(D) = 0, then the 
exotic D c has the same baryon and lepton number as 
the Standard Model d c field. Now R P (D) = +1 and 



R P (D) = — 1 as with the quarks of the Standard Model. 
In this case both (|2.3p and (|2 .4[) are forbidden, leaving 
only the first two lines of (|2.2p . At the renormalizable 
level, therefore, this is an accidentally conserved quantum 
"D-number" for the exotic fields and they are stable. Op- 
erators connecting the fields D and D c to the Standard 
Model may be allowed at the non-renormalizable level, 
however, depending on the U(l)' charge assignments. 6 In 
particular, for the case of the U(1)n combination, where 
B E = tan" 1 \/l5 0, E3, IH3], the combination of B, L, 
and U(1)n symmetry forbids the renormalizable opera- 
tors beyond the first two lines of (|2.2p . but allows the 
dimension-five operators 

dim5 : D c QH d S, D c QQu c , D c QLv c , (2.6) 

which preserve R-parity. These, along with the term pro- 
portional to A 5 in (12. 2p (which leads to a D,D C mass), 
allow for the decay of the exotics D and D c , which are 
therefore quasi-stable. An alternative model of quasi- 
stable exotics, in which a U(l)' gauge symmetry alone 
forbids D decay at the renormalizable level, can be found 
in Appendix [A] 

Finally, the case of mixing between the exotic D and 
SM d-quark leads to decays such as D — » uW, D — > dZ, 
and in some cases to D — > dZ' , or D — > d+ Higgs, where 
the W, Z, Z', or Higgs can be real or virtual. Such mix- 
ing can be induced by the operator \ s Dd c v c in the pres- 
ence of a sneutrino vev. Such examples are often consid- 
ered in the case of extensions of the MSSM in which one 
assigns L(v c ) = 0, as is often put forward in rank-6 mod- 
els. Mixing can also be induced by the operator X 7 D C QL 
if the scalar component of the neutrino in L acquires a 
vev, or by other operators such as D c QH ( i that are not in- 
cluded in (|2.2p because they don't occur in the singlet of 
27 3 and therefore violate the extra U (1) symmetries. The 
case of mixing between exotics and SM quarks (through 
arbitrary mechanisms) and its phenomen olog y has been 
well-covered in the literature [53, [5J, [55|, |j56| , especially 
in the non-supersymmetric case. We will therefore focus 
on cases where such a sneutrino vev or other mechanism 
is absent for the rest of this work. 



C. Mass Patterns 

We list in (12. 2p the operators A 4 and A 5 for complete- 
ness, but they are not particularly relevant for our study. 
The field S is a singlet under the Standard Model gauge 
group, but generally carries charges under additional 
Z7(l)'s which arise from the breaking of Eq to the Stan- 
dard Model. A vev for this field would generate both an 
effective fi parameter as well as a supersymmetric mass 



Fundamental bilinears are not allowed by Eq gauge invariance 
with only fundamental 27 representations. Furthermore, if we 
imagine a string-theoretic origin for our exotic D and D c states 
then such terms are generally forbidden if these fields are to be 
considered part of the massless spectrum of the string. 



6 The possibility of strongly interacting or charged exotics that are 
absolutely stable, or which are stable on the time scale of the age 
of the universe, was commented on in the Introduction. 
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for the D and D c . It will also generally break one or more 
additional t/(l)'s, producing new (heavy) Z'-bosons. As 
these facts are not relevant to the phenomenology we 
will pursue in subsequent sections, we will not consider 
them further. We do note, however, that the scalar mass 
matrices for the D and D c will generally depend on the 
charges of these fields under any additional U(l)' through 
various D-terms. 

Like the squarks and sleptons of the MSSM, the exotic 
scalar sector can also be described by a 6 x 6 scalar mass 
matrix. Unlike the MSSM fields, however, the fermionic 
modes are not protected by Standard Model chiral sym- 
metries from receiving large masses. In fact, we generally 
expect the scalar and fermionic components to receive a 
common, supersymmetric mass determined by the vev of 
some field such as the singlets Si in (|2.2p . Let us for the 
moment restrict our attention to the case of one gener- 
ation of exotics. Defining (S) = s and keeping in mind 
the definitions A 4 s = fi e g and A 5 s = van we can write 
the scalar mass matrix as 



111; 




(2.7) 



with 



in;, 



m 2 ~ + m 2 D + - sin 2 9 W cos 2/3 M% 



D 
J 2 



g' 2 Q' D (Q' s s 2 + Q 



H u V u 



YhA) 



>bi> 



m~ c + m D skr 9 W cos 2/3M| 

3 

+9' 2 Q' D AQ' s s 2 + Q'hA + Q'hA) 



where ( H' 



(2.8) 



= v Ut d are the usual scalar Higgs vevs 

and tan/3 = v u /vd- The quantity A§ is the soft 
supersymmetry-breaking A-term associated with the 
Yukawa interaction A 5 in (|2.2[) . We have ignored pos- 
sible CP-violating phases. In this work, equations (|2.7p 
and (|2.8[) are the only places where the precise Q 1 charges 
of the fields will be required. Note that the D-term con- 
tributions to m\ a and m\ h are typically small perturba- 
tions when g' ~ gy, where gy = y/5/3gi is the GUT- 
normalized U (l)y coupling, and can be absorbed into the 

values of the soft masses to- and m~ at low energies. 

d D c ° 

The resulting masses for the U(l) v model and U(1)n 
model are given for five benchmark points in Table |Tl] 
We will use these benchmark points to illustrate aspects 
of collider phenomenology in Sections IIVI and |V] below. 

For a particular set of U(l)' charges we can define the 
mass splittings 



A 2 



m Dl/2 
m Dl/2 



m D i 



l D 2 



(2.9) 
(2.10) 



between the physical fermion -D1/2 (with mass to_d 1/2 = 
mo) and the lightest scalar Dq or heavier scalar Dq, 



Parameter 


A 


B 


C 


D 


E 


M Dl/2 


300 


300 


300 


600 


1000 


m Do 


400 


400 


1000 


400 


400 


rriDg 


400 


400 


1000 


400 


400 


A B 


350 


150 


100 


600 


1050 



U(l) v Model 





367 


441 


1024 


388 


318 




587 


553 


1053 


932 


1482 


U(1)n Model 


M D i 


360 


435 


1022 


381 


309 


M D 2 


582 


527 


1050 


929 


1480 



TABLE II: Sample spectra for the exotic SUSY sector. 

Five benchmark mass patterns designed to illustrate the pos- 
sible collider signatures of exotic supermultiplets. All values 
are in GeV at the electroweak scale. These examples will be 
used extensively in Sections IIVI and IVl below. 



respectively. These splittings are functions of the dimen- 
sionful parameters mo, mjj, to^ c , A§ and /i e ffi as well as 
the dimensionless parameters g' , tan/3 and A 4 (or alter- 
natively the ratio v u Vd/s 2 ). If the quantity Ai is positive, 
then the lightest exotic particle (LEP) is the scalar. If it 
is negative then the LEP is the fermion. 7 As an exam- 
ple of the types of mass hierarchies that are possible in 
this parameter space, let us take the charges of the U(l) v 
model with g' = gy at the electroweak scale. Let us also 
fix the value of tan/3 = 10, // e a = 370 GeV and A 4 = 0.5 
(implying s = 740 GeV). 

Contours of Ai and A 2 for fixed values of the (com- 
mon) scalar mass mrj = = m^ a = 400 GeV are 
shown in Figures [T] and respectively. The shaded re- 
gions in the lower right of the plots are theoretically ex- 
cluded in that they produce a tachyonic eigenvalue of the 
exotic scalar mass matrix. For these low values of the 
soft scalar masses the majority of the parameter space 
admits a hierarchy in which the scalar is the LEP. In 
the limit as to-d — > (as is effectively the case for the 
fermions of the Standard Model), or in the limit where 
Ac, — > 0, the hierarchy is such that the scalars are gener- 
ally heavier than the fermions. This region has a lower 
bound dictated by direct searches limits on the mass of 
exotic SU(3) charged fermions, which we will discuss in 
the next section. For convenience we have labeled the 
relevant benchmark points from Table HT1 

In Figures [3] and U the same pair of quantities are plot- 
ted but with the supersymmetric fermion mass mo held 
fixed at 300 GeV. Again, the region in the lower right is 
excluded from theoretical grounds. From these plots it is 
clear that the fermion is the LEP unless the soft scalar 



7 Note that with the conventions of 112.101 1 the quantity A2 will 
generally be negative. 
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FIG. 1: Mass hierarchy Ai for fixed soft scalar masses. The 

quantity Ai = WlDj/j — m r> 1 * s plotted as a function of the fermion 
mass rri£) and the soft supersymmetry-breaking trilinear A5 for 
fixed (common) soft scalar masses mo = = mg c = 400 GeV. 
The shaded region in the lower right produces a negative mass- 
squared for the lightest scalar exotic. Relevant points from Table ITT1 
are labeled. 
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mion masses. The quantity Ai 



1/2 



m D i is plotted as 
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D. Production at Hadron Colliders 



masses are smaller than (or on the order of) the fermion 
mass and/or the trilinear j4-term coupling is sufficiently 
large. 



We expect strongly-interacting particles with the 
masses given in Table [TT] to be produced relatively fre- 
quently at hadron colliders. In this section we will dis- 
cuss the production cross-section for the leptoquark and 
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diquark cases before considering the direct search limits 
in the next section. Some aspects of the production of 
exotic SU (3)-charged states have been considered else- 
where [H, [H, [I?], [H, IH, [6(| HH, at varying levels of 
sophistication and approximation. 

As our goal is to be as complete as possible, we 
will consider the following ten 2^2 production pro- 
cesses: qq -f D 1/2 D 1/ 2, gg -> D 1/2 D 1/2 , qq -> D D , 
gg — ► D D and qg —> D + / (and c.c), with five 
each for the leptoquark and diquark cases. In addition, 
the couplings of Wjjq in (|2.4[) allow for the possibility 
of resonant production of exotic diquark scalars through 
quark or anti-quark annihilation. Where unavailable in 
the literature (or where available expressions were in- 
complete) we have computed the relevant parton-level 
production cross-sections to leading order and checked 
them against the results of CompHEP [62| . These expres- 
sions have been collected in Appendix B. The numeri- 
cal evaluation of these cross-sections - as well as all col- 
lider analysis performed in this work - was carried out 
with the PYTHIA 6.327 computer package [HI]- While 
the publicly-available version of PYTHIA does contain a 
scalar leptoquark, it does not have its superpartner, nor 
the diquark cases we wish to study. In addition, the 
scalar leptoquark contained in PYTHIA does not interact 
with the fields of the MSSM and can only decay into a 
quark and a charged lepton. Therefore some substantial 
modification to the off-the-shelf PYTHIA package was re- 
quired. We wish to briefly describe these modifications 
here in this section before proceeding. Further details of 
the analysis tools will be given in Section fVl 

Adding the desired new particles and interactions re- 
quired the modification of three existing subroutines and 
the addition of three new routines. Six new particles (two 
scalars and a fermion for the leptoquark and the diquark) 
were added to empty positions in the relevant common 
blocks, specifically the PYDAT2, PYDAT3 and PYDAT4 com- 
mon blocks. Masses and mixings of the new states were 
computed using the formulae of (|2.8[) via a new routine 
which parallels that of PYTHRG for standard MSSM sca- 
lars. A call to this new routine was inserted into the 
pre-existing PYMSIN SUSY initialization subroutine. De- 
cay rates for the exotics into Standard Model and MSSM 
states are computed and the necessary decay tables popu- 
lated with a new subroutine which is called from PYINIT. 
We will discuss the specific decay products considered in 
Section [V] below. 

The eleven new production processes were inserted into 
empty positions in the relevant common blocks, namely 
PYINT2, PYINT4 and PYINT6. The parton-level cross- 
sections were computed in a new subroutine called from 
the PYSIGH master routine. The most significant modi- 
fication of a pre-existing routine involved PYSCAT, which 
sets up the hard scattering process and documents the 
color flow through the interaction. For the leptoquark 
interaction, standard PYTHIA color flow algorithms suf- 
fice, but not so for the diquark interactions of (|2.4[) . 
These vertices involve three triplets or three anti-triplets 
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FIG. 5: Production cross section for pairs of leptoquarks 

at the LHC. Pair production of exotic fermions (g g,qq —* D D c ) 
is given by the solid (red) contours, while that of scalars is given 
by the dotted (black) contours. The region of coupling A ^ 0.2 
suggested by the indirect constraints considered in Section II II I is 
indicated by the light shading. 



of SU(3) - an interaction not present in the Standard 
Model. Such cases were not part of the original menu 
of color flow options in PYTHIA, so new ones were de- 
signed and inserted into the IC0L array for both diquark 
pair production and resonant production of scalar di- 
quarks. The essence of these modifications was to gener- 
ate place-holding "junctions" to serve as sinks or sources 
of color/anti-color. This modification is in the spirit of 
those used to study R-parity or baryon-number violating 
processes in the MSSM [11]. 

The above modifications allow us to simulate the 
eleven hard-scattering processes at LHC energies. For 
the sake of simplicity we will always take A 6 = A 7 and 
A 9 = A 10 in performing simulation-based calculations. 
We will refer to this common coupling as A, understand- 
ing that a different A value is implied for the leptoquark 
and the diquark. Resonant production of scalar diquarks 
was studied in detail elsewhere [65L [66| ; we postpone dis- 
cussion of this case to Section [Vj The production cross- 
sections for leptoquarks are given in Figures [5] and [6l 
while those for the diquark case are given in Figures 
and [5J Pair production cross-sections of exotic quarks 
and squarks are given in Figures [5] and [7| as a function 
of the mass of the exotic particle (denoted collectively as 
Mjy) and the Yukawa coupling A. Exotic scalar produc- 
tion in association with a Standard Model fermion via the 
process qg — > Do / is shown in Figures [JJ] and [51 In all 
figures we have shaded the region of small Yukawa cou- 
pling A < 0.2. In Section Hn] we will see that this may be 
taken as a very crude estimate of the allowed values of a 
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FIG. 6: Production cross section for scalar leptoquarks 
in association with fermions at the LHC. Contours give the 
production cross section for the process qg — » Do f. The region of 
coupling A 0.2 suggested by the indirect constraints considered 
in Section IlIII is indicated by the light shading. 



FIG. 8: Production cross section for scalar diquarks in 
association with fermions at the LHC. Same as Figure|6]but 
for diquarks. 
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essary to specify the masses of the superpartners of the 
Standard Model fields in order to unambiguously com- 
pute the production rate at the LHC. For the analysis 
presented here we will choose the well-studied bench- 
mark model SPS 1A from the "Snowmass Points and 
Slopes" collection [67|], in which the relevant masses are 
m Si ~ m fil = 535 GeV and m Sl = 146 GeV. The full 
set of superpartner masses for this benchmark point will 
be discussed in Section IVl below. 

The rate for production of exotic fermions is roughly 
an order of magnitude larger than that for identical-mass 
scalars, as one typically expects |3J. The five cases in Ta- 
ble |TT] were specifically chosen to give at least one exotic 
state in the 300-400 GeV range, ensuring a reasonable 
production rate at the LHC. In fact, the total production 
rate of exotics in Cases A-C in Table |H] is roughly equiv- 
alent to the total production rate of "standard" MSSM 
superpartners for the SPS 1A benchmark model. This 
implies that it should be possible to place meaningful lim- 
its on exotic masses and couplings from direct searches 
at existing colliders. We therefore turn our attention to 
direct and indirect experimental constraints on these pa- 
rameters. 



typical Yukawa coupling in this class of models. As these 
bounds are sensitive to many model-dependent phenom- 
ena we have chosen to display the cross-sections over a 
wide range of Yukawa parameters. 

Pair production of exotic fermions via the process 
q q — > D1/2 D1/2 can proceed through ^-channel exchange 
of scalar quarks and/or scalar leptons. It is therefore nec- 



III. CURRENT EXPERIMENTAL BOUNDS 
A. Direct Search Constraints 

The exotic quarks D and D c are charged under SU (3) 
and (as we demonstrated in the previous section) can 
thus be produced in large numbers through QCD pro- 
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cesses. Limits can be placed on their masses by mea- 
suring the cross-section x branching ratio for the exotic 
scalars and fermions into certain final-state topologies. 
These branching fractions - and to some extent the pro- 
duction rates as well - depend on the values of the allowed 
Yukawa interactions such as those in (|2.2[) . Here we wish 
to briefly summarize the limits placed on certain mani- 
festations of exotic SU(3) triplets from various collider 
searches. 

For diquarks current limits extend only to scalars 
which decay exclusively to pairs of jets. These jets can 
be produced resonantly with an exotic scalar in the s- 
channel. The CDF experiment was able to exclude di- 
quarks of the Eq type at the 95% confidence level in the 
mass range of roughly 300 to 450 GeV by measuring the 
cross-section to produce a resonant pair of dijets in a 
certain invariant mass window [H, EH, Iz3] • For reduced 
branching ratios (which is the case for much of our pa- 
rameter space) the limit essentially disappears. 

Leptoquarks have been more thoroughly studied at a 
number of collider environments. The DO experiment 
reported a limit in Run I of the Tevatron on the pair 
production of scalar leptoquarks which decay to the fi- 
nal states wqq, Ivqq and itqq for first and second gen- 
eration charged leptons. This corresponds to a mass 
limit at the 95% confidence level of 98 GeV for ex- 
clusive decay to the quark/neutrino final state. For 
Br (Do — > lq) = 1 the reported limit was rn D i 200 GeV 
and for Br(Do — > £q) — 0.2 a limit of m D i 150 GeV was 
given (jj, [ill, [73| . These results were updated at Run II 
and combined with the Run I data. The Run II results 
give m D i > 136 GeV at 95% confidence level for pure 
v q decays [7^ . The combined limits at the 95% confi- 
dence level give m D i > 256 GeV for Br(£>o — *■ eq) = 1 
and m D x > 234 GeV for Br(L> -> eq) = 0.5 [z|. 
Similarly, for second generation couplings the limits are 
m D i > 251 GeV for Br(Z?o — > fiq) = 1 and m D i > 
204°GeV for Br(L> M) = 0.5 [76]. DO also searched 
for third generation scalar leptoquarks, obtaining [73] 
m D i > 229 GeV for Br(D -> ^b) = 1. A simi- 
lar analysis was performed at CDF in Run II for sca- 
lar leptoquarks decaying to ilqq and ivqq final states. 
For first generation couplings the corresponding limits 
are m D i > 236 GeV, 205 GeV, 145 GeV, 126 GeV 
for Br(I>o -► eq) = 1, 0.5, 0.1, 0.01, respectively [z|. 
For second generation couplings the limits are m D i > 
226 GeV, 208 GeV, 143 GeV, 125 GeV for the same 
branching fractions to charged muons [7^]. Finally, DO 
searched for the production of a second generation lep- 
toquark in association with a fi. Combining with their 
associated production results, they obtained [801 ] m D i > 
274 GeV for A = 1 and Br(L> -> w) = l - 

Scalar leptoquarks can be produced on resonance 
at HERA. Here the production rates depend on the 
strength of the interactions in (|2.3|) . as do the rel- 
ative branching fraction to charged and neutral lep- 
tons. Assuming equal branching fractions, the limit from 



the ZEUS experiment is m D i > 290 GeV for A 9 = A 10 = 
A = 0.3 and m D x £ 270 GeV for A = 0.1 [Uj]. The limits 
from HI are similar (82| . 

Quasi-stable exotics require more specific search 
strategies which will depend on the manner in which 
they interact with the elements of the detector. The best 
limits come from the CDF search for massive charged 
hadrons which interact only weakly with the calorimeter 
but are tracked in the muon system. The lack of observed 
events puts a limit on the production cross-section for 
such exotic hadrons which corresponds at the 95% con- 
fidence level to an exotic fermion of charge \q\ — 1/3 of 
mn <; 190 GeV [HI]. Much weaker bounds on hadrons 
built from squarks and gluinos have been obtained from 
ALEPH at LEP 



B. Indirect Bounds 

There are a great many constraints on R-parity vi- 
olating operators in the MSSM 8 . In addition to direct 
searches at colliders, there are stringent indirect con- 
straints from proton decay (which essentially forbid the 
simultaneous presence of diquark and leptoquark opera- 
tors), neutron oscillations, K — K and B — B mixing, CP 
violation, rare B decays, lepton number and lepton fla- 
vor violating processes, neutrino masses, cosmology and 
astrophysics, and many other sources. 

Many of these processes also constrain the leptoquark 
and diquark couplings A 6 — A 10 defined in (|2.3[) and (|2.4|) 
of the exotic supermultiplets D and D c . As described 
in Section IIIB1 we assume that either the leptoquark 
operators or the diquark operators are present, but not 
both, and also that the scalar neutrinos v and v c do not 
acquire vacuum expectation values. In that case there 
are conserved baryon and lepton numbers and R-parity, 
implying the absence of proton decay (at least from the 
terms in (|2.2p ) and neutron oscillations, and also that 
there is no mixing between D or D c and the ordinary d 
or d c quarks. There are nevertheless many constraints 
from rare processes, analogous to those in the MSSM 
with R-parity violation, involving an internal D and/or 
D c line. These were studied many years ago by Campbell 
et al. [HI assuming specific relations between the scalar 
and fermion exotic masses and the other supersymme- 
try breaking parameters. A reanalysis of the constraints 
with more recent experimental values and general mass 
parameters is beyond the scope of this paper. 9 Moreover, 
there are many couplings involved when family indices 
arc included, often allowing individual constraints to be 
evaded by judiciously tuned choices of the dominant ones, 
and there is also the possibility of cancelations between 



8 For recent reviews, see [85l. |86|| . Older reviews include l87l |88|. 

Him mi. 

9 Some specific processes have been considered in [93, 94]. 
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diagrams. Because of these uncertainties, we will simply 
utilize the most stringent MSSM constraints for which 
there is a simple correspondence with a relevant diagram 
for exotic particle exchange, with the understanding that 
the constraints are to be considered a rough guide rather 
than a rigorous limit. Furthermore, our focus is mainly 
on couplings to the first two generations. Weaker bounds 
typically apply to couplings to the third generation. 10 

The most stringent relevant constraint on the lepto- 
quark couplings is from the SINDRUM II limit 



a(fx~Ti - 



e-Ti) 



capture) 



< 4.3 x 10" 13 



(3.11) 



on /i — e conversion [95( . There are several relevant tree- 
level diagrams in the R-parity violating MSSM, but only 
s-channel Dq or Dq exchange are relevant in the exotic 
model. From the MSSM analyses [H, [§7[ we estimate 



I \6 \6* 



< 



8 x 10~ 8 ( 



KH) Ciev) 



or 



|A 6 | < 3 x 10~ 4 ( 



m Da 



100 GeV 



(3.12) 



(3.13) 



if we assume \\% u c 



A 6 . A similar con- 



straint applies to A with m p —> rape . 

The limits on diquark couplings arc much weaker. The 
most important indirect limit involving the first two gen- 
erations (and ignoring possible CP-violating phases) is 
from the Kl — K$ mass difference. Th is has been an - 
alyzed in detail for the MSSM [H, M, DM EqU, IT02T ]. 
for which there are important contribution from box dia- 
grams involving four diquark vertices and also from boxes 
involving a W exchange as w ell as two diquark vertices. 
Using the estimates of |l00j . the most important dia- 
grams for the exotic case involve boxes with two internal 
exotic supermultiplet lines and two internal u, c, or t 
supermultiplet lines. Again ignoring the possibilities of 
cancelations or fine-tuning the family indices, one finds 
that typically 



A 9 ' 10 < 0.04 



max(m Si ,m,D J \ 



100 GeV 



(3.14) 



for the couplings involving external d or s and internal 
u,,. 



C. Cosmological Bounds 

A number of significant constraints on the proper- 
ties of quasi-stable particles, particularly quasi-stable 



10 For that reason the analysis in [2tJ assumed that the leptoquark 
or diquark couplings only involved the third generation. 



hadrons, arise from cosmological observations. The most 
severe bounds come from the successful prediction of 
light element abundances by the Big Bang Nucleosyn- 
thesis (BBN) model. When the value of the baryon-to- 
photon ratio r\ = ub/tIj of (6. 10 ± 0.21) x 10~ 10 from 
the WMAP three year results 103] is used, the theo- 
retical predictions of the standard BBN scenario are in 
reasonable agreement with the observed abundances of 
deuterium D, 3 He, 4 He, 6 Li and 7 Li. These successful 
predictions are based, however, on a set of assumptions 
about early-universe physics. In particular they assume 
the Standard Model field content and set of interactions. 
The presence of late-decaying exotic particles beyond the 
Standard Model are therefore likely to change the physics 
that ultimately gives rise to the primordial abundances 
of light elements. 

Constraints on long-lived exotics arising from BBN de- 
pend on two key quantities: the abundance of the exotic 
and its lifetime. Shorter lifetimes are generally less well- 
constrained than longer ones, but even if the exotic par- 
ticle decays wells before the epoch of BBN (beginning 
roughly at the temperature scale T ~ 0.7 MeV), there 
can be conflict with the successful BBN predictions if 
the abundance of the exotic is too large. The relic abun- 
dance is a function of the annihilation cross-section for 
the (un-hadronized) exotics in the early universe. We 
will consider the issue of relic abundance below, but first 
we summarize the principal observational constraints on 
any new long-lived particles. 

If some new state decays with a lifetime greater than 
roughly a tenth of a second, then it will decay during 
or after the epoch at which BBN takes place, and the 
decay products can potentially alter the predictions of 
BBN. They can mediate additional interconversion be- 
tween protons and neutrons beyond that of the Standard 
Model interactions at early stages of the BBN process. At 
later stages these decay products can cause hadrodissoci- 
ation or photodissociation of the primordial background 
nuclei. The implications of these n on-standard processes 
are explained in detail in jl04l . Il05l | . We will here merely 
sketch the conclusions before applying these results to 
the model of Section [TTJ 

Constraints on new long-lived particles apply to states 
with lifetimes r in the range 10~ 2 sec ^ t ^ 10 12 sec. 
Roughly speaking, the type of constraint depends on the 
epoch in which the new state decays, and this can be 
divided into three temporal regions. For states decay- 
ing with lifetime 10~ 2 sec r ^ 10 2 sec hadronic de- 
cay products from the exotic are likely to lose energy 
very quickly through electromagnetic processes. They 
are therefore insufficiently energetic to destroy the newly- 
created light element nuclei. Nevertheless, strong interac- 
tions allow scattering of these decay products with pro- 
tons and neutrons, generating interconversion between 
the two (in particular conversion of protons into neu- 
trons). This is even after these baryons freeze out with 
respect to electroweak conversion processes. The result- 
ing increase in the ratio of neutrons to protons results 
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FIG. 9: BBN Constraints on New Hadronically Decaying 
P artic les. This figure, reprinted with permission from the authors 
of [105| , summarizes the constraints arising from BBN on a particle 
of mass mx = 1 TeV that decays exclusively to hadronic Standard 
Model states. The various contours represent bounds arising from 
the observation of different primordial element abundances. The 
vertical axis indicates the relic density Yx (i.e., the ratio of number 
density to entropy) of the exotic multiplied by the amount of energy 
released into interacting particles per exotic decay. For this plot 
that fraction is 100%, or E v i s = mx- 



in a higher yield of deuterium and 4 He than is observed 
experimentally. 

For longer-lived particles with 10 2 sec ^ r ^ 10 7 sec 
the mesons from the decay products tend to decay be- 
fore they have a chance to interact with the background 
neutrons and protons. Thus the n/p ratio is likely to 
be unchanged. But the hadronic decay products from 
the exotics will be emitted with a higher kinetic energy 
than the thermalized nuclei of the light elements. In this 
case hadrodissociation of alpha-nuclei is common and the 
result is nonthermal production of deuterium and 6 Li. 
Finally, for the case of very long-lived exotics in which 
10 7 sec t ^ 10 12 sec even the neutrons from the decay 
of the exotics will now have time to decay before inter- 
acting with background nuclei. In this case photodisso- 
ciation from emitted photons and hadrodissociation are 
competitive processes. The constraint arises from non- 
thermal overproduction of 3 He. 

These constraints from BBN are summarized in Fig- 
ure [H where the region above and to the right of the var- 
ious curves are excluded. 11 As indicated in the preceding 
paragraphs, the bounds arise from different observations 
depending on the lifetime r of the exotic. For shorter 



lifetimes the principal bound arises ultimately from ob- 
servations of primordial 4 He (as indicated by the mass 
fraction Y p ) and primordial deuterium (as indicated by 
the ratio D/H). The pairs of curves for these two ob- 
servations reflect different estimates for extracting the 
primordial abundances from current observations. For 
the case of the primordial mass fraction o f 4 H e these es- 
timates are from Izoto v an d Thuan (IT) [1061 ] and from 
Fields and Olive (FO) [l07j |. In the case of deuterium the 
"high" and "low" estimat es di ffer in whether or not the 
recent data of Webb et al. [l08j | is included in the fit. For 
intermediate lifetimes the 6 Li observations are most con- 
straining, while very long lifetimes are most constrained 
by observations of 3 He. 

The relic abundance is determined by the annihilation 
cross-section. For definiteness we will assume that the 
lightest exotic particle (LEP) is a fermion. 12 The an- 
nihilation proceeds most often through QCD processes 
into quarks and gluons, with a thermally-averaged cross 
section 



{a\v\) 



(14 + 6N f )ir 
27 



D 



(3.15) 



where the first term is for gluons and Nf is the number of 
quark species that may appear in the final state. For our 
purposes we will take Nf — 5. This, and the correspond- 
ing formula for scalar LEPs, 13 can easily be computed 
using appropriate modifications of the cross sections for 
the inverse processes listed in Appendix [B] In addition, 
we will assume that 50% of the energy is carried away 
by non-interacting decay products (such as the LSP of 
a supersymmctric theory). With these assumptions, the 
nature of the BBN constraint will depend crucially on 
the exotic lifetime. 

Let us consider the model suggested in section III Bl in 
which all renormalizable operators allowing for the de- 
cay of the exotic SU(3) triplets/antitriplets are forbid- 
den. Then the first case of allowed decay operators arise 
at mass dimension five via the operators listed in (|2.6p . 
There are therefore three kinds of decay channels for the 
quasi-stable fermion D, depending on its mass. If the 
D is heavier than the two superpartner bosonic particles 
and/or the sum of the scalar Higgs and S masses, then 
three body decays into a fermion and two scalars will 



It is also quite possible that the mass difference between the sca- 
lar and fermion is smaller than the LSP mass. In this instance 
both the fermion and scalar would be quasi-stable. Similar state- 
ments apply to heavier generations of exotics. 
If the LEP is a scalar, the averaged annihilation cross section 

(2»+6N f T/m DQ )TT 



at temperature T is (<t|^|) 



where 



This plot is reprinted from Ref. |105|l with permission of K. Kohri. 



the T/mo facto r is because the annihilation into quarks is p- 
wave suppressed |l09j . Since T/m£> ~ 1/20 at freezeout, the 
annihilation is mainly into gluons, with a cross section about 2/3 
that for the case of a fermion LSP, and the corresponding relic 
density about 50% higher. 
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dominate. 14 Should this decay be kinematically forbid- 
den, two-body decays into a massless fermion and scalar 
Higgs or scalar S may be induced through the first opera- 
tor of (|2.6p . As we expect our exotic states to have masses 
well in excess of 100 GeV to avoid the direct search con- 
straints of Section IIII A[ we will assume that the decay 
channel D — > dh is always available through this oper- 
ator. Lastly, the same operator will allow the decays 
D — ► dZ, uW through mixing effects. If the exotic fer- 
mion D is lighter than all scalar quarks and leptons it 
will decay primarily through these processes. 

For three body decays into a (massless) fermion and 
two (massive) scalars the decay width is given by 



C, 



5127T 3 



x K* 



(3.16) 



where Ci is a numerical factor that depends on the par- 
ticular final state and K3 represents the integral over 
phase space. It is 1/3 neglecting the scalar masses, and 
in general is 



(2 — xa — XB)dxA<ixi 



(3.17) 



where x a.b — 2EA t B/mD with Ea,b the energies of sca- 
lars A and B. The integration limits are given in the 

Appendix of [lid] . The quantity M* is the mass scale 
that suppresses the higher-dimensional operators in the 
superpotential. Two body decays into a fermion and sca- 
lar are given by 



Ci 



D 



(h d ) 



64tt 



m 2 D Ml 



x K 



s,v 



where 



= 1 



""D 

2mi 
rrin 



1 - 



D 



(3.18) 

(3.19) 
(3.20) 



for decays into scalars (jua = mh or m s ) and vectors 
(iriA — mw or mz), respectively. 

Clearly the resulting lifetime t& of the exotic fermion 
is a strong function of the relative masses of the states in 
the theory and will be sensitive to precisely how many de- 
cay channels are available to the exotic. To be concrete, 
let us consider the 17(1) jy case described in Section HlBl 
with the dimension-five operators of (|2.6|) . Only the first 
operator allows for mixing or two-body decays (assum- 
ing no vevs for the superpartners of the Standard Model 
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FIG. 10: Lifetime tjj of the Exotic Quark in the Quasi- 
Stable Scenario of Section IIII The lifetime as a function of 
exotic mass is plotted for various values of the suppression scale 
AT. of 10 19 , 10 18 , 10 17 , 10 16 , 10 15 and 10 14 GeV. Abrupt changes 
in the curves are the result of additional decay channels opening as 
the mass Mjj increases. 



fields or v c ). A conservative estimate of the lifetime can 
be made by assuming only decays into first generation 
fields, allowing only those channels for which the fermion 
final states are massless (i.e., no decays into Higgsinos 
or singlinos), and assuming no decays into heavy Higgs 
states or Z'. This results in 14 partial widths for the fi- 
nal states: dZ, uW, dh, ds, dsh, dLU^u, dLU R u, dv,LU R , 
u L e L v, u L ev R , ue L v R , d L v L v, d L vv R , dv L v R . The pref- 
actors for these cases are given by 



Cdz = 



4G F 



Cdh = 



V2 ' 
4GV sin 2 a 
V2 cos 2 /? 



C u w = 



Cds = 



8G F 
1 



Cdsh = sin 2 a/4 ; C q qq = 4 , 



(3.21) 



and unity for the final states arising from the final oper- 
ator of (|2.6p . a is the usual MSSM mixing angle for the 
Higgs scalar mass eigenstates. For the decays involving 
the second operator in (|2.6[) . with only first generation 
particles, the two fields arising from the two quark dou- 
blets (of which there are two possible combinations) must 
not have the same color. Hence the total counting factor 
gives Ci — 4. 

For completeness, we also consider the decays medi- 
ated by virtual W's and Z's for mr> ^ my/ (decays in- 
volving virtual h, s, or Z' are much less important than 
other allowed decays for all masses). The result of a 
straightforward calculation, using the partial rates cal- 
culated in [531, is 



14 The smaller phase space for three body decays is compensated 
by the larger number of channels, by color factors, and by the 
cos 2 /3 suppression of the two body rates, at least for the specific 
model considered here. 
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is the D L — d L mixing angle (assumed small) and 
sin 2 8w ~ 0.23 is the weak angle. In (|3.22p we have as- 
sumed that the D mixes only with the d, neglected CKM 
mixing, and included decays into three families of mass- 
less leptons and five flavors of massless quarks. W — Z 
interference and identical particle effects are included for 
the duu and ddd channels, respectively. 

To obtain numerical estimates it is necessary to postu- 
late specific mass values for certain scalar fields. Let us 
take rrij = 500 GeV for all scalars of the MSSM, m s = 
1000 GeV for the singlet scalar field, 15 m h = 115 GeV, 
(s) = 740 GeV and tan (3 = 10. We also will take the de- 
coupling limit for the MSSM Higgs sector such that the 
mixing angle a is given by a — f3 — it/2. The resulting 
lifetime is plotted in Figure ITD1 as a function of exotic fer- 
mion mass Mo and the mass scale of the dimension five 
operator M* . As the actual numerical values of these re- 
sults are sensitive to the masses of the scalars (as well as 
the number of decay channels considered), the results in 
Figure [TTj] should be taken as indicative of what is likely 
in models of this sort. For values of Mo < rrih decays 
can only proceed via mixing with Standard Model fer- 
mions; thus the lifetime decreases slowly with increasing 
Mr>, Most of this range can generally be excluded from 
direct collider searches (though this is a model-dependent 
statement) as discussed in Section ITlI Al Above this mass 
scale the lifetime generally decreases rapidly, particularly 
for Mr, greater than twice the typical scalar mass of the 
squarks and sleptons. We conclude from Figure [T0l that 
for reasonable exotic masses and suppression scales M* 
(i.e., M* less than or equal to the reduced Planck mass) 
the primary constraints will involve proton-neutron con- 
version and hadrodissociation of background alpha nu- 
clei. 

From the cross-section (|3.15l) the freeze-out tempera- 
ture Tf and relic density Yd of the exo tic p article can 
be computed using standard techniques [109j . For each 
value of M* there exists a unique pair of values for to 
and E vm Yo (or to and Mo) which can be compared with 
the experimental bounds. In Figure [TT] we have plotted 
these contours as a function of M* for the values 10 18 , 
10 17 , 10 16 , 10 15 and 10 14 GeV and overlaid them on the 
limits from Figure [TJ] The darker shaded region is the 
union of all constraints from Figure [9] using the weak- 
est bounds for D/H and Y p , while the lighter shaded 
region extends this disallowed space by using the more 
restrictive values for these quantities. Along each curve 



15 We will only consider the real part of this scalar field to be dy- 
namical. When the models of this paper are embedded in theories 
with additional C/(l)'s it is the imaginary part of this field that 
is typically "eaten" to produce a massive Z' boson. 
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FIG. 11: BBN Constraints on Exotic Quark Parameters. 

We plot combinations of lifetime rr> and Mrj (or equivalently 
E via Y D for NU of 10 18 , 10 17 , 10 16 , 10 15 and 10 14 GeV. The darker 
shaded region is the union of all constraints from Figurc[9]using the 
weakest bounds for D/H and Y p , while the lighter shaded region 
extends this disallowed space by using the more restrictive values 
for these quantities. Note that this plot assumes that 50% of the 
mass energy of the exotic is carried away by non-interacting LSPs 
(such that E vis = 0.5M D ). 



the value of Mo increases from very small to very large 
values as one moves from the lower right to the upper 
left. For large suppression factors (M* ^ 10 17 GeV) only 
extremely light exotics (or extremely massive ones) can 
be tolerated by existing limits. Such small values would 
likely be in conflict with collider bounds from the Teva- 
tron. For M* ;$ 10 15 GeV nearly the entire range of mass 
values Mo is allowed by current observations. 

We again stress that there are several factors that can 
change the impact of the light element abundance con- 
straints on the exotic mass. Some of these factors would 
serve to weaken the bounds on Mo for a given M*: 
changing the scalar masses for superpartners (in partic- 
ular m s ), allowing additional decay channels into Higgsi- 
nos, singlinos and second/third generation quarks, allow- 
ing additional energy to be carried off by non-interacting 
particles, and so forth. Furthermore, it has recently been 
suggested that hadronized exotics should be expected to 
undergo a second round of annihilation shortly after the 
QCD phase transition 111]. This can serve to dilute 
these bounds considerably as well. We therefore con- 
clude that it is not unreasonable to assume that should 
string-inspired exotic states exist, with decay amplitudes 
suppressed by a mass scale somewhat below the reduced 
Planck scale, then they may very well exist in mass ranges 
relevant to the upcoming LHC experiment. Such scales 
could emerge in string constructions with some dimen- 
sions larger than the Planck length, or be associated with 
non-Planck scale physics. In the next section we will con- 
sider the phenomenological implications of this scenario, 
before returning to the case of prompt decays in Sec- 
tion El 
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IV. QUASI-STABLE EXOTICS AT THE LHC 

Given the discussion from the previous section, it is 
reasonable to assume that exotic new triplet/anti-triplet 
SU(3) representations will hadronize into color-singlet 
states on a rapid time scale. These exotic hadrons have 
come to be known as R-hadrons in the literature and we 
will adopt that name for them here. These R-hadrons 
will be stable on the time scale of the detector and can 
be treated as hadrons for the duration of their interaction 
with the detector elements. There has been some discus- 
sion of such quasi-stable exotics in the literature, but this 
has been primar ily w ith r egard to lo ng-li ved gluino-based 
R-hadrons [HI Ell El LLlS EH, E3 ■ While the phe- 
nomenology of long-lived triplet/anti-triplet representa- 
tions will be similar, there are important differences. We 
will therefore revisit the subject in this section and focus 
on the signature of such states at the LHC experiments. 
Our analysis will be sufficient to give a sense of the dis- 
covery reach of the LHC for such new states, though a 
more refined analysis with a full detector simulation will 
no doubt sharpen the conclusions obtained here. 



A. Collider Phenomenology 

Unlike the previous studies which assume a gluino com- 
ponent for the R-hadron, here the exotic component can 
be a scalar or a fermion. Furthermore, these earlier 
papers were generally motivated by scenarios with very 
heavy scalars, thus assuming only direct pair production 
of gluinos. Here there is the very real possibility of pro- 
ducing pairs of the heavier states in the supermultiplet; 
the next-to-lightest exotic particles (NLEPs). Depending 
on the mass differences involved, these states may then 
decay into the LEP which subsequently hadronizes. This 
provides an additional handle for triggering and event 
reconstruction which we will describe in the next subsec- 
tion. 

Another important distinction between the case of the 
gluino and that of the exotic quark is in their SU(3) rep- 
resentation, which impacts the sorts of R-hadrons that 
can form and their interaction cross-sections with the 
detector elements. Like their gluino counterparts, ex- 
otic triplets and anti-triplets can form R-mesons and R- 
baryons, with the former kinematically favored at the ini- 
tial hadronization stage. However, the R-hadrons consid- 
ered here will necessarily have one fewer "active" quark 
than those of the gluino-based variety. As a result, the 
total cross-section for R-hadron interactions with protons 
and neutrons will be proportionally reduced. 

The R-meson states will include Dd and Du combi- 
nations (as well as the anti-states) and the R-baryons 
will include Ddd, Duu and two combinations of Ddu 



(as well as their anti-states). 16 The meson states will 
be approximately degenerate in mass, with the bulk 
of the mass being accounted for by the exotic compo- 
nent. We thus expect roughly 50% of the R-mesons 
formed at the initial vertex to be charged. Once pro- 
duced, charged R-mesons will leave tracks in the inner 
detector elements. The R-mesons will pass largely unaf- 
fected through the electromagnetic calorimeter and enter 
the hadronic calorimeter. Here the R-hadrons will un- 
dergo a series of interactions which include elastic scat- 
tering off nucleons, charge-exchange interactions with 
nucle ons a nd meson-to-baryon/baryon-to-meson interac- 
tions [114]. These processes yield an interaction cross- 
section of roughly 12 mb for R-mesons and 24 mb for 
R-baryons. The array of possible interactions depends 
on the R-hadron itself and reveals an important asym- 
metry. 

R-mesons of the form Dq will preferentially undergo 
meson-to-baryon transitions by extracting two quarks 
from a target nucleon and producing a light pion in the 
final state. The resulting Dqq R-baryon will remain an R- 
baryon due to the absence of anti-quarks in the detector 
material. Among the possible R-baryons which can form, 
the case Dud with d and u in an s-wave configuration will 
be the lightest state, with the p-wave configurations of 
Dud, Ddd and Duu more massive. In contrast, R-mesons 
of the form D c q will remain R-mesons due to the lack of 
quark-antiquark annihilation possibilities. Even if an R- 
baryon of the form D c qq were to form, it would quickly 
be destroyed by baryon-to-meson interactions which pro- 
duce light pions. 

With the cross-sections given above, we expect the ex- 
otic R-hadrons to undergo a number of interactions with 
the hadron calorimeter. For R-baryons (R-mesons) we 
estimate on average 8 (6) interactions in the "barrel" 
region (0 < n < 1.5) and 10 (7) interactions in the "end- 
cap" region (1.5 < n < 2.4) for both the ATLAS and 
CMS detectors prior to entering the muon system. As 
the bulk of the energy and momentum is carried by the 
(non-interacting) exotic quark, these interactions tend to 
resul t in a very small energy deposit in the calorimeter 
cells jll4| . The exact amount depends on the mass of the 
R-hadron, its kinetic energy and the material through 
which it is passing. As an example, consider the pair 
production via QCD of exotic fermion LEPs with mass 
Mo = 600 GeV. Figure [T2l gives the distribution of ki- 
netic energies of the produced LEPs from a simulation 
10,000 events (roughly 10fb _1 of data-taking). While 
most are produced with relatively little kinetic energy 
there is a long tail in the distribution , wit h a mean at 
E kin ~ 400 GeV. From the results of [115| we estimate 



Approximately 15% of the R-hadrons produced in the primary 
high-py process will involve a strange quark, with a resulting 
6 mb reduction in their interactions with protons and neutrons. 
As this will have a negligible impact on the collider physics con- 
sidered below we will neglect all heavy flavor R-hadrons. 
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FIG. 12: Initial Kinetic Energy of Produced R-hadrons. 

Distribution of kinetic energies for fermionic LEPs with mjj = 
600 GeV prior to entering the calorimeter system. 



the typical energy loss per interaction for R-hadrons with 
-Ekin 1 TeV to be somewhere between 0.2 and 2.2 GeV. 
For our simulations in the next subsection we will assume 
1.5 GeV in energy loss per nuclear interaction. 

We therefore expect most, but not all, R-mesons pro- 
duced in the primary interaction with \r}\ < 2.4 to 
punch-through to the muon chambers as some form of 
R-hadron in the ATLAS or CMS detector. Those includ- 
ing D will most often arrive as neutral R-baryons while 
those involving D c will arrive as some form of neutral or 
charged R-meson. If we neglect the small possibility of 
D c s R-mesons, we would anticipate roughly 50% of the 
R-mesons involving D c to be charged when they enter 
the muon chamber. These charged R-hadrons will leave 
tracks in the muon system and will have a characteris- 
tic velocity [3 significantly different from the essentially 
massless muon. Of course it is possible for R-hadrons 
of all types to further interact hadronically in the muon 
system - perhaps in a charge-exchanging way. We will 
neglect that small possibility in what follows. 

Typical distances between the widely separated resis- 
tive plate chambers in the muon systems of ATLAS and 
CMS are less than 1 meter, with a separation between 
the first and last such plate at approximately 3 me- 
ters |l!8j . It should therefore be possible to make a 
robust distinction between tracks arising from the ex- 
otics and those from background muons by measuring 
the respective times-of- flight (TOF). The primary source 
of background muons will be single and double weak bo- 
son production. Separating the signal from background 
can be performed by requiring that the TOF between 
any pair of reference points for the candidate R-hadron 
be at least 3 ns greater than the TOF of a muon (with 
/3 ~ 1). The value of 3 ns is sufficiently greater than the 
8t resolution of both the ATLAS and CMS m uon system 
to provide a highly significant S/y/B value [l!5| . The 
final requirement is that the R-hadron be moving suf- 
ficiently swiftly to arrive at the muon chamber with the 
time-window for the data to be recorded with the current 



B. Discovery Reach 

There remains, however, the question of triggering on 
these events. What fraction of the R-hadron events will 
eventually be recorded to tape? The low-level trigger 
must capture an event prior to full reconstruction, thus 
the presence of tracks in the inner detector will not be suf- 
ficient to capture an event without significant activity in 
either the calorimeters or the muon chambers. R-hadrons 
which stop in the calorimeter will likely deposit sufficient 
energy to trigger in the E^ m channel or i?™ lss channel 
(assuming the other R-hadron punches-through or decays 
at a different time). However, given the potentially very 
long lifetimes of these states these decays will occur at 
a much later bunch-crossing than that which produced 
the exotic quarks. Such issues have been addressed in 
the context of long-lived gluinos [120j . We will not pur- 
sue the phenomenology of these cases further, though we 
will note the number of such "stopping" exotics in our 
simulations to follow. 

The punch-through R-hadrons are likely to leave only 
10-50 GeV of transverse energy in the detector - not 
enough to trigger in the E^ m channel. Furthermore, as 
the exotic quarks are pair produced and are back-to-back 
in the center-of-mass frame, the total amount of trans- 
verse energy carried away by the exotics is small. Of 
course, production of the next-to-lightest exotic parti- 
cle (NLEP) allows for the possibility of either increased 
E™ lss through Dq — > -D1/2X1 decays or increased jet ac- 
tivity through Do — > D 1 / 2 g decays. In the former case, 
the two D1/2X1 systems are again back-to-back, with only 
a slight increase in E™ lss ; typical values are less than 
50 GeV. 

If either of these triggers is utilized, the event will al- 
most certainly pass the second level triggers as well (due 
to the presence of, say, a charged track leading to the 
calorimeter cells and/or the presence of a track in the 
muon chamber). Nevertheless, it is still preferable for 
a track to be identified in the muon chamber for signal 
extraction and correct particle identification. We thus 
consider the muon trigger alone. Since we expect the 
emergence of only one charged R-meson (at most) from 
the two produced exotics, the appropriate low-level trig- 
ger is the single-muon channel. We will require a very 
conservative pt threshold of 15 GeV for triggering on 
the R-mesons that enter the muon chamber (after ac- 
counting for the energy loss due to hadronic interactions 
in the calorimeter). 

We illustrate the effectiveness of this search strategy in 
Table IIIII for the five benchmark cases presented in Sec- 
tion UTTJ The acceptance at each stage is roughly con- 
stant across the benchmark scenarios. Approximately 75- 
85% of the R-hadrons are produced with \rj\ < 2.4, and 
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Benchmark Point 





A 


B 


C 


D 


E 


Geom. Accept. 


75.5% 


79.9% 


82.3% 


86.8% 


82.5% 


Charged Frac. 


25.2% 


25.0% 


25.1% 


25.2% 


25.4% 


Temp. Accept. 


82.7% 


82.8% 


81.9% 


79.1% 


76.9% 


TOF 


97.3% 


96.5% 


97.2% 


97.3% 


97.0% 


Total Accept. 


15.3% 


16.0% 


16.5% 


16.9% 


15.6% 


iVsignal (Xl0 3 ) 


120 


119 


119 


11.2 


26.6 


iVstop (xlO 3 ) 


11.1 


10.8 


11.3 


1.36 


4.56 



TABLE III: Signal Acceptance for Quasi-Stable R- 
hadron Scenarios. Geometrical acceptance represents the 
fraction of R-hadrons that are produced with |?7| < 2.4. 
Temporal acceptance represents the fraction of charged non- 
stopping R-hadrons that arrive within 18 ns of the primary 
interaction for the event. The percentage that traverse a 3 me- 
ter fiducial distance at least 3 ns slower than a /3 = 1 muon 
would is given by TOF. The product of these fractions is the 
total acceptance. The number of signal events (as well as 
the number of stopping R-hadrons) is given for 10 fb _1 of 
integrated luminosity. 



roughly 25% emerge from the calorimeter into the muon 
system as charged mesons. Of these, approximately 80% 
have p > 0.5 and thus arrive within 18 ns of the pri- 
mary interaction in the event. Each of these R-hadrons 
will therefore produce a charged track in the muon sys- 
tem. The distribution in transverse momentum for these 
objects upon arrival at the muon chambers is given in 
Figure IT51 for Scenario C. In this case all of the R- mesons 
have sufficient pt to trigger given our 15 GeV minimum 
Pt requirement. This fact was true of all five benchmark 
points. Thus adding additional trigger possibilities (such 
as E™ lsa ) is unlikely to add significant numbers of sig- 
nal events if the muon system is to be used for particle 
identification. Finally, the fraction of R-mesons moving 
sufficiently slowly to traverse a 3 meter fiducial distance 
at least 3 ns longer than a f3 = 1 muon is given by the 
"TOF" entry in Table HTT1 This represents the vast ma- 
jority of R-mesons that enter the muon system within 
the 18 ns time window. We therefore estimate the total 
acceptance to be approximately one-sixth of all produced 
quasi-stable exotics. 

The discovery reach will track the production cross- 
section for the lightest exotic particle. In each of our 
benchmark cases there is at least one exotic state with a 
mass below 450 GeV, providing for copious production at 
LHC energies. For the simulations described in Table Hill 
we allow for production of both the exotic fermion and 
the lighter exotic scalar. In case A both the fermion LEP 
and its slightly heavier scalar will be quasi-stable, while 
the other cases will involve SUSY cascade decays for the 
heavier states. The higher production cross-section for 
the fermion (as demonstrated in Section Hi Dp reflects it- 
self in the factor of ten between the number of signal 
events that arise for cases A-C and those of cases D and E. 
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FIG. 13: Transverse Momentum of Charged R-mesons. 

The distribution of px for charged R-Mesons with /3 > 0.5 upon en- 
tering the muon system. We assume a minimum of px > 15 GeV 
to trigger on the charged track. All R-hadrons moving with the 
minimum velocity have sufficient momentum to meet this thresh- 
old. 
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FIG. 14: Reach in Exotic Mass for Muon Signature. Dis- 
covery reach at the ATLAS experiment for "punch-through" quasi- 
stable exotics. The threshold for discovery is taken to be ten muon- 
like events, assuming a negligible background rate. 



Given the exceptionally large signal-to-background ratio 
for events of this type, discovery will not prove a problem 
if such light exotics exist. We also note that a fair number 
of the produced exotics have insufficient kinetic energy to 
punch-through to the muon system and will stop in the 
calorimeter. The issu e of detecting these events has been 
addressed elsewhere [lll| . As the lifetime of these ex- 
otics is an undetermined parameter, we merely list the 
number of such events in Table IIIII 

As the masses of the exotics increase, the production 
rate falls. We estimate the discovery reach in nifl and 
TO£) via the muon signature channel in Figure 1141 We 
use the standard practice of taking ten events to consti- 
tute the discovery threshold in cases with an exception- 
ally small number of expected background events. This 
corresponds to a mass reach of tod ^ 1700 GeV for the 
fermion and tod ^ 1450 GeV for the scalar exotic in 10 
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FIG. 15: Stopping R-hadrons as a Function of Exotic Mass. 

Discovery reach at the ATLAS experiment for quasi-stable exotics 
which stop in the calorimeter and subsequently decay. The thresh- 
old for discovery is taken to be ten such events, assuming a negli- 
gible background rate. 



fb -1 of integrated luminosity. Since it should also be able 
to detect the fraction of events that stop somewhere in 
the detector itself, we also include the mass reach for this 
signature in Figure [151 Here the mass reach for 10 events 
is 771 u 1100 GeV for the fermion and tod 925 GeV 
for the scalar. These rough estimates for the discovery 
reach of quasi-stable ST/(3)-charged objects are in good 
agreeme nt with calcula tions performed in other model 
contexts [HI [HI El]. 



V. PROMPT DECAY SIGNATURES 

Our intention is to study the presence of exotic super- 
multiplcts in the presence of "standard" supersymmet- 
ric states such as those of the MSSM. This is motivated 
by a number of considerations. First, the study of the 
collider signatures of _E@-inspired exotics in isolation has 
been performed by a number of authors. The signatures 
and analysis strategies for these cases are by now well 
understood. For example, a scalar state which decays 
directly into a lepton + quark or into two quarks can be 
reliably reconstructed by looking at the invariant mass 
distribution of leptons + jets and pairs of jets, respec- 
tively. This is possible even in the presence of Standard 
Model backgrounds and a measurement of the mass of 
the exotic scalar can readily be made (66| . It is therefore 
more interesting to consider the situation when a richer 
spectrum of "new physics" is present. 

Second, most theoretical models of supersymmetric 
physics that have an ultraviolet completion - whether 
a GUT model or some sort of string construction - tend 
to predict that the exotics described in Section |TT] will 
be present in conjunction with the states of the MSSM. 
In fact, as alluded to in the Introduction, they will typi- 
cally arise with additional gauge fields (and their gaugino 
partners) and additional scalars. Therefore, considering 
the exotics in isolation may oversimplify the experimental 



challenges associated with these fields. Finally, it is inter- 
esting to ask how the presence of MSSM superpartners 
can serve as "background" to the identification, extrac- 
tion and eventual study of exotic new supermultiplets - 
and indeed how the exotic states can complicate other- 
wise well-established analysis routines for supersymmet- 
ric parameter measurements. 

In this section we will study the collider signatures 
of our exotic benchmark cases from Table [TTJ We will 
focus primarily on the leptoquark scenario for technical 
reasons, but we will occasionally address the diquark sce- 
nario as well. We reserve a full treatment of the diquark 
cases to a future publication. We allow the full set of 
interactions between exotics and MSSM states given in 
equations (|2.3p or (|2.4p , but restrict ourselves to interac- 
tions between the exotic states and the first generation 
of Standard Model/MSSM states. We will also restrict 
the interactions with gauginos and Higgsinos to include 
only those states of the MSSM (i.e., we continue to ne- 
glect the possible presence of additional singlets or U(l ) 
factors). This truncation reduces the model-dependence 
of the study considerably and corre sponds to plausible 
limits of full-scale extended models [l2ll Il22l | . We will 
also restrict attention to a single exotic fermion and the 
lighter exotic scalar, and will ignore the possibility of 
heavier generations of exotic particles. The latter would 
generally be expected to decay rapidly to the lightest gen- 
eration if they mix and are sufficiently heavy, or else to 
exhibit production and decay patterns similar to to those 
of the lightest generation. Finally, we ignore decays into 
the right handed neutrino v c and its scalar partner. De- 
pending on the model of neutrino mass generation, their 
masses might or might not be relevant at the LHC [4(| . 
If they are sufficiently light, the effects of the A 8 term in 
(|2.3[) would be similar to those from the neutrino part of 
the A 7 term. 



A. Decays of Exotics 

The majority of previous studies on the collider phe- 
nomenology of leptoquarks and diquarks have focused 
on non-supersymmetric scenarios in which the exotic 
state decays directly to Standard Model states. This 
may arise through mixing operators or by consider- 
ing only the production and decay of scalar exotics 
with branching fractions to pairs of SM fermions set to 
unity [H, [56|, [U [H, [H, [6(| [6l| . In order to discuss the 
full range of supersymmetric decays available - and to 
study the interesting issues outlined in the introduction 
to this section, it is necessary to postulate a fixed spec- 
trum of MSSM states to include with the five scenarios 
of Table [Hi 

Our MSSM "background" sample corresponds to the 
minimal supergravity point SPS la, from the Snowmass 
Points & Slopes J67j. This point has a spectrum deter- 
mined from the following set of minimal supergravity pa- 
rameters: a unified scalar mass ttiq = 100 GeV, a unified 
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Par am. 


SPS la 


Par am. 


SPS la 


m~o 
m~o 


99.9 
188.4 


m h 

™>h 


381.4 
587.3 


rn~a 

*3 


375.5 




535.3 


m~o 
x 4 

m~± 
x± 

m~± 


394.0 
187.7 
394.7 
627.9 


m-r 

771s, . 777,7 


554.5 
504.5 
535.0 
534.4 


B-ino% 


97.4% 




559.3 


m h 


111.7 




145.5 


rn A 


412.7 




220.6 


mf 


420.3 




145.8 


A* 


369.4 




211.4 



TABLE IV: Sample non-exotic spectra for Snowmass 
point la (SPS la). All masses are in GeV. 
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FIG. 16: Relative mass scales for SPS la and exotic cases 
of Table lilt The mass difference between the exotic fermion and 
the lightest scalar, given by the quantity Ai in 1)2. 9jl is given for 
each of the five cases summarized in Table [TT] We indicate with 
dotted arrows the masses of some of the relevant MSSM states 
that may appear in the decay chains. 



gaugino mass TO1/2 = 250 GeV, a unified trilinear scalar 
coupling A Q = —100 GeV, tan/3 = 10 and positive fi pa- 
rameter. These soft supersymmetry-breaking values are 
specified at a high-energy input scale, here assumed to 
be the GUT scale (A GUT = 1 x 10 16 GeV). These GUT- 
scale parameters are passed directly to PYTHIA, which 
then generates the physical mass spectrum for the MSSM 
at the electroweak scale. This process involves using 
approximate analytical solutions to the renormalization 
group equations (RGEs) for the soft Lagrangian param- 
eters. Superpartner masses are computed at tree level, 
except for the Higgs sector which is computed from the 
full one-loop effective potential. The result is given in 
Table ITVl A more careful treatment of the RGEs and/or 
the physical eigenstate masses might correct these num- 
bers by as much as 10%, but the accuracy achieved here 
is sufficient for out purposes. 

Snowmass Point la is the benchmark point from (6?| 
most favorable for SUSY discovery at the LHC. It is 
therefore one of the most-studied supersymmetric models 
in the literature. All the superpartners are within reach 
of the LHC (and indeed, some may be detectable at the 
Tevatron with sufficient integrated luminosity) . The rel- 
atively light scalar leptons are particularly important for 
this benchmark point, as they significantly enhance the 
event rate for trileptons and large missing transverse en- 
ergy ($t) through the production of electroweak gaugi- 
nos. The relatively light gluino also ensures a large rate 
for multi-jet events with large We will discuss these 
properties in the following subsections. 

The five cases in Table [TT] were selected to provide a 
range of possible decay chains between the exotic states, 
the MSSM spectrum and Standard Model states. The 
most important of the relevant mass scales are shown 
schematically in Figure [TH] Mass differences between the 
exotic fermion and the lightest exotic scalar are given in 
terms of the Ai variable introduced in (|2.9[) for the five 



cases (note that Ai < for cases A-C). Key thresholds 
are indicated by superpartner masses for the two lightest 
neutralinos x? 2> the gluino g, the lightest first-generation 

slcptons e± and the lightest first-generation squarks d±. 
Note that the exotic fermion and its lightest scalar part- 
ner are each above 300 GeV in mass in all five cases. 

In the first three cases of Table HTJ the fermion is the 
lightest exotic particle (LEP). For the leptoquark sce- 
nario the two-body decay channels to a scalar lepton 
and a (massless) quark are always available due to the 
very low-mass sleptons in the SPS la model. Taking 
A 6 = A 7 = 0.1, the branching fractions for into 

the final states d v% , it and u are 28%, 50% and 
22%, respectively for cases A-C. For the diquark scenario 
the equivalent two-body decays to a quark and squark 
arc kincmatically forbidden. In this case the fermionic 
diquark must decay via three-body channels into states 
involving gauginos and two SM quarks. For our analysis 
we will only consider cases in which the final state in- 
volves the LSP neutralino. These three-body decay rates 
are given in Appendix B for decays via first-generation 
virtual squarks, as are those for the leptoquark case (rel- 
evant in different kinematic regimes for which the two- 
body decays are forbidden). 

For both the leptoquark and diquark scenarios the di- 
rect decay into Standard Model fermions is the only chan- 
nel kinematically available to the exotic scalar in Case A. 
The diquark has a single channel — > u d, while the 
leptoquark has two final states — ► u e~, dv e with 
relative branching fractions of 67% and 33%, respectively. 
For Case B the channel Dq — ► Di/ 2 Xi opens up, but the 
decays to fermion pairs continue to dominate. In the 
leptoquark case Br(Z?o Q -> D\^ 2 x?) = 15% with the 
two SM channels sharing the remainder in the ratio of 
2:1. For the diquark case Br(D° Q -► x?) = 5% for 
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FIG. 17: Branching ratios into various final states for Dq® 
decays in Scenario C. The branching fraction into the seven 
allowed final states for the scalar leptoquark are given as a function 
of a (universal) Yukawa coupling A 6 = A 7 = Alq. These allowed 
final states are as follows: (1) D 1 / 2 g (2) -D1/2X? (3) ue (4) dv e 
(5) D 1/2 xl (6) D 1/2 xi (7) D 1/2 xt 



FIG. 18: Branching ratios into various final states for -D^/a 
decays in Scenario D. The branching fraction into the eight 
allowed final states for the fermionic leptoquark are given as a 
function of a (universal) Yukawa coupling A 6 = A 7 = Alq- These 
allowed final states are as follows: (1) Do xi (2,3,4) du e , ueu, u&l 
(5) eu R (6) eu L (7) ud L (8) D x°2- 



A 9 = A 10 = 0.1, with the other 95% accounted for by the 
two quark final state. For Case C all four neutralino final 
states are available, as is the decay process Dq — > D 1 / 2 <?■ 
This latter mode is the dominant channel, with a branch- 
ing fraction of 87% in the leptoquark case and 77% in the 
diquark case. The next largest branching fraction for the 
leptoquark is Br(L>g Q -> D\9 2 x?) = 8% and for the di- 
quark is Br(Dg Q — > ud) = 16%. The decay possibilities 
for the cases A-C are summarized in Table fVl As an ex- 
ample we give in Figure [17] the branching fraction into 
the seven allowed final states for decays of the scalar D]^ 
as a function of the size of the leptoquark Yukawa inter- 
action. For convenience we here take A 6 = A 7 = Alq. 
The vertical dotted line at Alq = 0.1 defines the model 
we will consider in Section IV Bl 
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partner + xi 
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partner + X3 












/ 


partner + xi 
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partner + g 












// 


f + f 
f + f 


//lq 
NA 


NA 
// 


//lq 
NA 


NA 
// 


//lq 
NA 


NA 
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Xi + .f + f 


//dq 




//dq 




//dq 





TABLE V: Decay possibilities for cases A-C in Table JT] A 

checkmark indicates an allowed decay channel. The double check- 
mark indicates the channel with the largest branching fraction. 
Note that this channel need not be the same for the diquark and 
leptoquark scenarios. 

The scalar partner is the LEP for cases D and E. Here 
the scalar always decays via the Yukawa couplings of (|2.3j) 
and (|2.4p to pairs of SM fermions. Case D was designed 



to allow the decays Di/ 2 —> D X121 as weu as decays via 
the Yukawa interaction to / /' final states. The branch- 
ing fraction for decay into the exotic scalar + LSP is 32% 
for the leptoquark and 72% for the diquark. In the latter 
case, decays into qq' final states are suppressed by the 
relatively large squark mass scale. Final states such as 
d v e and u e~ are similarly suppressed for the leptoquark. 
In this latter case the final states dv^, u e ^ and u e£ 
each have branching fractions in the 20-24% range, mak- 
ing them competitive with the Dq \i final state. The 
branching fraction for D^® 2 into the eight allowed final 
states are plotted in Figure [18] as a function of Yukawa 
coupling Alq- Again we indicate the choice Alq = 0.1 by 
the vertical dotted line. For Case E the very large mass 
difference between the fermion and scalar LEP allows all 
decay processes for both leptoquark and diquark. Here 
the dominant rate is to Dq g with branching fractions of 
97% and 94%, respectively. The decay possibilities for 
cases D and E are summarized in Table [VI] 
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/ 




partner + xX 
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f + f 
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/ 

NA 


NA 
// 


/ 

NA 


NA 
// 



TABLE VI: Decay possibilities for cases D and E in Ta- 
ble [ij A checkmark indicates an allowed decay channel. The 
double checkmark indicates the channel with the largest branching 
fraction. 

Each individual decay channel was given its own po- 
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sition in unused entries in the BRAT and KFDP common 
blocks in PYTHIA. Two-body decays are computed at ini- 
tialization by the new routine PYEXDC on the basis of 
user-designated Yukawa interactions and particle masses. 
Three-body decays require a more sophisticated treat- 
ment due to the numerical phase space integrations re- 
quired. For these specific decays the result of these nu- 
merical integrations are given directly to PYTHIA for the 
five scenario points in Tabic [TTJ 



B. Collider Signatures 

We now turn our attention to the signature for these 
exotic supermultiplets at the CERN LHC. We arc inter- 
ested in what can be learned about the exotic component 
of the SUSY signal in the earliest stages of LHC data col- 
lection and analysis. Hence we will study 5 fb _1 of sim- 
ulated signal data, where the signal here includes both 
the events associated with the SPS la benchmark point 
and one of our five exotic scenarios. Signatures involv- 
ing the diquark tend to be multi-jet events with slightly 
less missing transverse energy than a "typical" MSSM 
SUSY event. As we will not be simulating Standard 
Model backgrounds for QCD jet events, making meaning- 
ful statements about signal-to-background estimations is 
difficult, even with a full detector simulation. A realistic 
study of such exotics would require a more sophisticated 
treatment of SM multi-jet backgrounds than we are pro- 
viding in the current work. Furthermore, an additional 
technical difficulty arises in this case when hadronizing 
the out-going partons in PYTHIA when non-standard (i.e. 
non-planar) color flows are involved. 17 For these two rea- 
sons we will defer the treatment of the case of diquark 
couplings for the exotic multiplets to a subsequent publi- 
cation, when the technical issues mentioned here can be 
more adequately resolved. 

For the leptoquark scenario the Standard Model back- 
grounds are similar to those of typical SUSY events and 
have been studied elsewhere jl23l . 1124 , Il25l | . In this pa- 
per we will rely on those estimates to allow us to discuss 
the collider signatures without computing the relevant 
background samples. Nevertheless, a sample of 200,000 
tt and 40,000 W + W~ /W ± Z/Z Z events were gener- 
ated for examining the efficiency of our cuts in isolat- 
ing the exotic leptoquark signal from both the Standard 
Model and MSSM backgrounds. We proceed by simulat- 
ing the e vents with PYTHIA and then pass the events to 
PGS [l26l | to simulate the detector response. The detector 
parameters are those of the hy brid ATLAS /CMS detec- 
tor used by the LHC Olympics [l27j | . As we do not wish 
to concern ourselves unduly with detailed issues of trig- 



Unfortunately, we have yet to find an adequate solution to this 
problem, which (to the authors' knowledge) afflicts all publicly- 
available event generation software. 



Sample 


Accepted Events 


Percentage 


oro la, 


1 8^ ^A A 


no 1 0/ 

yo. i /o 


Case A 


161,284 


85.0% 


Case B 


156,020 


84.7% 


Case C 


152,342 


84.5% 


Case D 


11,589 


96.8% 


Case E 


17.921 


94.0% 



TABLE VII: Number of events in 5 fb _1 of integrated lu- 
minosity. For our baseline model SPS la and each of the five 
leptoquark cases 5 fb — 1 of data was simulated at the LHC. The 
second column gives the number of these simulated events for which 
the H^ lg variable was over 150 GeV. The third column gives the 
fraction of generated events which pass this crude trigger threshold. 

gering, we have chosen to utilize the simplified triggering 
algorithm from version 3.0 of the PGS package. More 
specifically, we will analyze events for which the quantity 
i4 rig > 150 GeV, where 

^5xJ2 b£ P |+0.2x£ \p^\+0.2xJ2 b J T Ct M Vt\ ■ 

lep tau jet 

(5.24) 

The objects appearing in the summation above are can- 
didate objects; e.g. the electrons and muons in the first 
term in (|5.24j) are trigger-level electrons and muons. The 
taus in (|5.24|) are hadronically-decaying tau-candidates. 
Crudely speaking, electron and muon candidates must 
have pt > 10 GeV, hadronically-decaying tau candidates 
and jet candidates must havepy > 100 GeV and_^T must 
be greater than 50 GeV to appear in the trigger function. 

1. Can the new physics be definitively detected? 

As mentioned in the previous subsection, the spectrum 
for Snowmass Point la gives rise to signatures that have 
been extensively studied in the past. The presence of ex- 
cesses above Standard Model predictions can be defini- 
tively established in multiple channels in this scenario 
for even a handful of inverse femptobarns of integrated 
luminosity. The most important of these channels are 
multi-jet final states, multi-jets plus one or more leptons, 
events with same-sign pairs of leptons and any number 
of jets, trilept qn events , and events with three taus in 
the final state [l28l Il29l Il30j | . All of these cases are for 
events accompanied by large $t values. 

We expect the additional contribution from the iso- 
singlet exotics to add to these event rates - making the 
discovery of "new physics" even more rapid. In 5 fb _1 of 
integrated luminosity 185,544 MSSM events from SPS la 
passed our low-level trigger function (an acceptance rate 
of about 93%). The corresponding number of accepted 
events for the five exotic scenarios is given in Table IVIII 
The relatively large production rate for the light exotic 
fermion in cases A-C nearly doubles the new physics 
sample, while the light-scalar scenarios provide only a 
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FIG. 19: Effective mass distribution for multi-jet (non- 
leptonic) events. The effective mass variabfc (1 5 - 2 5 P > for the four 
hardest jets in events with no leptons and^y > 100 GeV. The five 
curves represent SPS la together with each of the five leptoquark 
cases. The bulk of the events passing the selection criteria are from 
SPS la, resulting in almost no distinction between SPS la and the 
five exotic scenarios. 
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FIG. 20: Effective mass distribution for all high-pr ob- 
jects in events with jets plus leptons. The same variable as 
Figure Il9l but now formed from all high-py objects with at least 
two jets and any number of leptons. In this case the additional 
contribution from pair production of exotic fermions is apparent in 
fermion-LEP cases A-C. 



small additional contribution. Not all of these additional 
events, however, will enter the SUSY search analyses. 

Let us consider the bluntest, but most effective, tool for 
establishing the presence of beyond-the-Standard Model 
physics: the M c g analy sis fo r events with mul tiple high- 
Pt jets and large ,Et .123]. Following [l25j ] we select 
events with at least four jets, no isolated leptons with 
Pt > 20 GeV, transverse sphericity S > 0.2, and fix > 
100 GeV. We require the pr of the hardest jet to be at 
least 100 GeV and the pr of the four hardest jets to each 
be above 50 GeV. The variable M c g is given by the scalar 
sum 

M eS = p Tjl + p Ta + Pt,3 + Pta+ #t, (5.25) 

where pT,i is the transverse momentum of the i-th jet, 
ranked in order of px value. We plot this quantity for the 
SPS la model together with each of the five exotic sce- 
narios in Figure [TO1 For mass bins with M e g <; 600 GeV 
a clear excess over the rapidly-falling Standard Model 
background should be visible. 

Note, however, that the additional exotic component 
adds very little to the event sample in this analysis. 
While 33,836 events in the SPS la sample pass the above 
event selection cuts, less than 3,000 events pass them in 
cases A-C, while for cases D and E only 410 and 490 
events, respectively, enter the analysis. Even if we allow 
events with any number of high-p-p leptons we increase 
the SPS la acceptance by a factor of 1.5, but for the 
exotic samples the increase is only a modest factor of 
3 for the fermion LEP cases and 2 for the scalar LEP 
cases. The most severe cuts turn out to be the require- 
ment of four hard jets and the cut on fix- We prefer 
not to relax the $t > 100 GeV cut without a thorough 



understanding of the Standard Model backgrounds. But 
in Figure [20] we plot the quantity M e g for all objects in 
events with at least two jets, each with pj, > 50 GeV 
and any number of leptons (the other cuts remaining the 
same as before). The scalar LEP cases continue to con- 
tribute very little to the analysis, but now some structure 
is evident in the M e s distribution for the fermionic LEP 
cases, albeit for low-energy bins where Standard Model 
processes will likely be dominant. In Figure |2"T1 we focus 
on Case A where we show the individual contributions to 
the M c ff variable of Figure [20] from SPS la and from the 
fermionic LEP Case A. The two distributions are clearly 
separated, with the SPS la peak representing squark and 
gluino production while the peak for exotic Case A rep- 
resents fermion LEP pair production. Whether or not 
the peak in this distribution at lower M e g values is ob- 
servable over the Standard Model background is an open 
question that will require further study with full back- 
ground simulation. 

As our exotic signal will generally involve one or more 
hard leptons in the final state we might expect that the 
inclusive signatures involving leptons will be more highly 
augmented over the SPS la case alone. In Figure [22] we 
give the ev ent r ate for th e six key inclusive signatures 
studied in [HI [TH, [Hj for SPS la and the five ex- 
otic scenarios of Table [n] The first inclusive signature 
involves the event sample for the M c g analysis: multi- 
jet events with large fir and no isolated leptons. The 
other five of these discovery modes involve high-py lep- 
tons. They are jet events with single isolated leptons, 
events with jets and opposite-sign (OS) dileptons, events 
with jets and same-sign (SS) dileptons, events with three 
leptons and jets and events with three high-p^ taus. All 
of these signatures involve large missing transverse en- 
ergy 
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FIG. 21: Breakdown of the effective mass distribution for 
exotic case A. The contribution of exotic states to the generalized 
M e ff variable is clearly separated from the contributions of standard 
MSSM squarks and gluinos. A more thorough analysis of Standard 
Model backgrounds is necessary, however, to conclude that this 
structure in the low-mass bins is experimentally observable. 
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FIG. 22: Number of events of different inclusive signa- 
tures. Six very inclusive signatures that may serve as "discovery 
modes" for supersymmetric new physics are presented. These sig- 
natures are: (A) multi-jet events, (B) jets + leptons, (C) OS dilep- 
tons, (D) SS dilcptons, (E) trileptons and (F) three-tau events. All 
signatures have at least two high-py jets and large $t- Note the 
logarithmic scale of the vertical axis. 



In contrast to Ref. Il30l we give numbers of events of 
these signature topologies after applying a set of cuts de- 
signed to reduce the Standard M odel back ground. These 
cuts are based on those of Refs. HH [HI All events in 
Figure [22] have fx > 100 GeV and n- Jct > 2, except for 
the multi-jet category where we require nj c t > 3. For the 
multi-jet case we require the three hardest jets to have 
Pj, > 100 GeV. For the other five signatures we require 
1?? > 100 GeV for the hardest jet and p£* > 50 GeV 
for the second-hardest jet. All events are required to 
have a transverse sphericity of S > 0.2. For the multi-jet 
case we veto events with leptons having \rj\ < 2.5 and 



Pt P ^ 20 GeV. For the topologies involving leptons we 
instead require that the leptons (or taus) be isolated and 
have p° v > 20 GeV. 

As expected, the multi-jet channel is overwhelmingly 
populated with events from the MSSM states. Event 
rates for the leptonic final states for cases A-C (with the 
fermion LEP) can be competitive with those of standard 
SUSY - particularly for the OS dilepton case where the 
large production cross-section for fermionic LEPs couples 
with the unit branching ratio for these states to quark + 
slepton. For cases D and E (with the scalar LEP) the 
exotic contribution to these inclusive signatures is a very 
small fraction of the total event rate. Nevertheless, in 
a world with only the supersymmetric exotic sector the 
scalar-LEP models are likely discoverable on their own in 
5 fb _1 of data in the trilepton and jets + lepton channels 
(particularly if judicious cuts on the size of $t ar e made 
in these events). The fermionic LEP models give strong 
signals in all leptonic channels. When combined with the 
large rate from SPS 1A it is clear that robust "discovery" 
modes for new physics will be present across the board, 
even in the early stages of LHC running. 



2. Can the exotic component be extracted? 

Observation of the combined event rates in the six 
channels of Figure [H] should give clear evidence of new 
physics at the LHC, and should furthermore strongly sug- 
gest supersymmetry as the source of that new physics sig- 
nal. Yet it is not immediately clear that the exotic com- 
ponent - above and beyond the spectrum of the MSSM - 
would be recognized merely from the inclusive signature 
counts. Taking into account the logarithmic scale in Fig- 
ure we note that the MSSM processes still dominate 
most of the inclusive signatures, particularly for exotic 
cases D and E. Nevertheless, the signal is clearly "lepton 
rich" relative to SPS la alone. Presumably an attempt 
to fit initial data to paradigms such as minimal super- 
gravity would be unable to reproduce the large rate for 
multi- lepton events. 18 The first suggestion of beyond- 
the-MSSM physics may therefore arise from fits to the 
data which suggest extremely light slepton masses (say, 
below the bounds set by direct search limits) or anoma- 
lously large branching fractions of gauginos to leptonic 
final states. 

One can find additional hints from the data itself. We 
have already seen that plotting the effective mass of all 
objects in events with at least 100 GeV of $t, at least two 
jets of pj? > 50 GeV and any number of leptons produces 
a structure indicating some subset of events are not aris- 
ing from gluino and/or squark production. Yet these low 



It would be an interesting exercise to attempt such a fit within 
the context of a suitably general MSSM framewo rk, s uch as the 
"minimal reasonable model" of Brhlik and Kane [ 1 3 lH ■ 
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FIG. 23: Missing transverse energy distribution. Total T^t 
in all events is plotted independently for the SPS la data sam- 
ple and the five exotic leptoquark cases. In fermionic LEP models 
(open histograms) the exotic fermions decay through cascade pro- 
cesses ending in the neutralino LSP. For scalar LEP models (filled 
histograms) there is very little T^t for scalar production, and very 
few events involving the much heavier exotic fermion. 



effective-mass bins may not be visibie above t he St an- 
dard Model background when the cuts of Ref. Il23l are 
loosened. Indeed, the rough topology of events arising 
from pair production of leptoquark D l / 2 pairs is not un- 
like top quark events; cuts designed to reduce the rather 
large tt background at the LHC will therefore tend to 
reduce the exotic signal in these channels as well. 

Nevertheless, some basic observations may suggest a 
leptoquark interpretation for at least some of the signal 
events. Those events arising from either fermionic or sca- 
lar leptoquarks distinguish themselves from those arising 
from the MSSM states of SPS la in three main ways: 
(1) lower jet multiplicities, (2) higher lepton multiplici- 
ties and (3) lower values of $t- In Figure [23] we plot the 
distribution of missing transverse energy in all events for 
each of our data samples. The standard MSSM processes 
give rise to a broad distribution which peaks at relatively 
large $t values. For example, requiring a $t cut of 100 
GeV only eliminates 9% of the data sample (as opposed 
to 75% of our t i sample and 87% of our diboson sample). 
In contrast, the fermionic LEP cases give a much sharper 
peak at roughly half the value of the SPS la case. Nearly 
50% of these events have fir < 100 GeV. For the sca- 
lar LEP cases we have a sizable production of the scalar 
exotic, which decays directly to two Standard Model fer- 
mions with almost no missing energy. The long shallow 
tails are from the production and decay of the much heav- 
ier exotic fermion. Missing transverse energy in Standard 
Model events tends to peak at or below 50 GeV and fall 
rapidly. This suggests that at least for the fermionic LEP 
cases, the bimodal distribution in J£t for the joint case 
of MSSM + exotics should indicate two sets of fields be- 
ing produced, both with production cross-sections sug- 
gestive of strongly-interacting states, but with different 
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FIG. 24: Lepton multiplicities. Numbers of events with one 
to five leptons, each having at least 20 GeV of transverse momen- 
tum, for SPS la and the five exotic leptoquark scenarios. Note the 
logarithmic vertical axis. 



mass scales and decay topologies. 

In Figure [23] we give the typical number of leptons per 
event for the individual data samples. In each case the 
leptons are required to have at least 20 GeV of transverse 
momentum, but no other cuts are applied to the events. 
The exotic leptoquark samples clearly favor higher lep- 
ton multiplicities, but the characteristics of these leptons 
themselves (such as M e g computed solely from leading 
leptons, or the pt distribution of the hardest lepton) 
do not differentiate between the SPS la events and the 
exotic events. The experimentalist investigating these 
events will find no peaks in the invariant mass distribu- 
tions of various pairs of leptons. However, the events 
with multiple high-pT leptons tend to be those events 
with less Pt- Consider, for example, exotic case A in 
conjunction with the SPS la model. In Figure [2FJ the 
scalar sum of p^ p values for all leptons in events with at 
least two isolated leptons is plotted versus the missing 
transverse energy in those events. It is clear from Fig- 
ure [55] that there is some rough anti-correlation between 
the energy accounted for by leptons and the amount of 
energy carried away by LSPs. The feature visible along 
the horizontal axis in the figure (that is, for extremely 
small values of $t) with extremely energetic leptons are 
events involving the pair production of scalar leptoquarks 
- which have a sizable event rate for the 367 GeV scalar 
mass of this benchmark point. The events in this part of 
the scatter plot also appear in the low fix bins in Fig- 
ure [23] for case A. All of the events in Figure [25] have 
high-pr jets, yet no such anti-correlation or low $t fea- 
ture is apparent in the analogous scatter plot involving 
jet pt values. 

At this stage, if we were confronted with cases A-C we 
might conclude that in addition to squarks and gluinos 
(with a mass-scale of approximately 600-800 GeV) we are 
also producing a second class of states with a mass-scale 
of approximately 300-500 GeV. The overall rate for pair 
production of new states of this mass is consistent with 
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| Case A: Events with two or more leptons I 
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FIG. 25: Scalar sum of lepton px values versus J9t for 
events in SPS la + exotic case A. Events in this data set 
were required to have two isolated leptons but no other cuts were 
made. The events at low values of Ifix and large lepton px values 
are events from exotic leptoquark production, particularly scalar 
pair production. 
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FIG. 26: Invariant mass of hardest lepton paired with soft- 
est jet in two jet, OS dilepton events. Precisely two jets, 
neither being B-tagged, were required, as were two opposite-sign 
leptons. For cases A C, with significant fermionic LEP produc- 
tion, the end-point in the jet-lepton invariant mass distribution at 
200 GeV correctly gives the mass difference Mp — Mlsp- For the 
four cases where scalar production was non-negligible a mass peak 
can be reconstructed near the physical mass value for the lightest 
scalar. 



those states being fermions which interact through the 
strong interaction. Events of this class occur with smaller 
amounts of fix than the squark/qluino events, and are 
associated with increased numbers of leptons and fewer 
jets. The association of strongly-coupled physics with 
prompt leptons and smaller but still significant J^t signals 
may suggest a supersymmetric leptoquark interpretation. 
For cases D and E, where the event rates are lower, the 
fix signal is muted or non-existent and the Standard 
Model backgrounds are more problematic; it is unlikely 
that any deviation from the MSSM spectrum would be 
apparent given the crude observations presented here. 

We have argued that 5 fb _1 of integrated luminosity is 
likely sufficient to reveal the presence of supersymmetric 
new physics. In the three cases involving a light exotic 
fermion it is also likely sufficient to reveal one or more 
objects with a production rate indicative of strong cou- 
pling but with decay chains in which high-py leptons are 
prominent. Even in cases D and E for which there will 
be little evidence for an exotic component to the SUSY 
signal, it is still likely that experimental searches for lep- 
toquarks will be performed nonetheless. 

In Figure we plot the invariant mass distribution of 
the hardest lepton and softest jet in events with precisely 
two jets and two (opposite-sign) leptons. Jets were re- 
quired to have at least 50 GeV of transverse momentum 
and events were vetoed if either jet was B-tagged. A cut 
was made on the pr of the leading lepton of 50 GeV, and 
20 GeV for the trailing lepton. Finally, we require the 
events to be somewhat collimated along the event axis, 
so we require the transverse sphericity to be no greater 
than 0.7. This final cut significantly reduced the contam- 
ination from both Standard Model processes and SPS la 
events (an acceptance rate of approximately 0.04% for 



each). 

For cases A-C the invariant mass of the jet/lepton pair 
shows an end-point just below 200 GeV. This correctly 
measures the mass combination 
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via the on-shell cascade decay -D1/2 — > ql —> qfyv F° r 
cases A-C this happens to be very near the mass differ- 
ence between the fermionic LEP and the lightest neu- 
tralino. Mass peaks arising from the scalar pair produc- 
tion with D — > q£ can be reconstructed for all scenarios 
in which there is significant scalar production (case C had 
only 38 scalar events in 5 fb _1 of data). The true mass 
value for the lighter scalar is given over the correspond- 
ing peak in Figure l26l We note that if a cut on missing 
energy of fir > 50 GeV were applied, the scalar mass 
peaks would vanish from the distributions in Figure |2"61 
though the end-point in the distribution associated with 
fermion pair-production would still be visible. 

These peaks can be isolated and sharpened by making 
stricter cuts on the data set, such as demanding fir < 
25 GeV, requiring the scalar sum of pt values from the 
two jets and two leptons sum to at least 400 GeV, and 
requiring the invariant mass of the lepton pair to be at 
least 100 GeV. An important cross-check is to find the 
same peak in the jet/lepton invariant mass distribution 
in associated production of scalar leptoquarks through 
the process gq — > D^q. We can isolate this process by 
requiring (a) at least two jets without B-tags, the hardest 
jet having at least 200 GeV of transverse momentum and 
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all others having px > 50 GeV, (b) precisely one isolated 
lepton with p T > 50 GeV, and (c) E T < 20 GeV. Pairing 
the second hardest jet with the single lepton gives a clear 
peak at the same mass values as those in Figure [5U 

Thus, in every one of the scenarios of Table [IT] there 
should be at least one exotic state, and occasionally two 
such states, which can be identified at the LHC - even 
with limited initial data. With additional statistics it 
should be possible to measure the masses of low-lying 
scalar mass eigenstates in all five scenarios. Reconstruc- 
tion of cascade decays with additional integrated lumi- 
nosity should also allow a determination of the exotic 
fermion mass in all five cases. We note, however, that 
the presence of the (supersymmetric) exotic states can 
complicate the analysis of the signal arising from MSSM 
states. 

As the simplest example, we may consider the flavor- 
subtracted dilepton invariant mass distribution for OS- 
dilepton events in Snowmass point SPS la. For low-mass 
SUSY points such as this, even a small data sample is 
sufficient to establish a sharp "edge" in the weighted sum 
of distributions given by 

M inv (e + e-) + M inv (fi+ f i-)-M inv (e + t i-)- M inv (fi+ e -) . 

(5.27) 

The quantity in (|5.27|) for SPS la is given by the dashed 
line in Figure [13 The edge at M inv (t+£-) ~ 90 GeV 
is the result of the flavor-subtraction procedure. When 
the pair of opposite-sign leptons arise from a single de- 
cay chain of x° — * t — > ^ + ^~X?, the two leptons are 
forced to be of the same flavor. When the two leptons 
come from two independent decays (as in the case of t — t 
production, or pair production of squarks and gluinos) 
the probabilities for same-flavor and opposite-flavor lep- 
ton pairs are roughly the same. The subtraction pro- 
cedure therefore efficiently isolates the exclusive decay 
chain for x§ decays and produces a sharp edge at the 
value 
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which happens to be very nearly the mass difference 
Mj — M^o in Snowmass point la. 

Yet the effectiveness of this subtraction procedure is a 
result of the flavor-democracy of Standard Model gauge 
interactions (which are related to the decay rates of gau- 
ginos by supersymmetry) . For our exotic states these de- 
cays are governed by Yukawa interactions which (up to 
model-dependent constraints from flavor-changing pro- 
cesses discussed in Section ITlI B[) are largely free parame- 
ters. When the fermion LEP preferentially decays to one 
flavor of lepton, for example, the subtraction technique 
will fail to completely isolate gaugino decays. We show 
this effect by the solid line in Figure |2"T1 where Snowmass 
point SPS la is analyzed together with the exotic states 
of Case A. In this case the exotic is decaying exclusively 
to first-generation leptons. The distribution still exhibits 
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FIG. 27: Flavor-subtracted dilepton invariant mass distri- 
bution. When flavor-universality is not present in the leptoquark 
Yukawa couplings flavor-subtraction will not be adequate to remove 
exotic fermion decays from x% decays. 



an edge (suggesting on-shell sleptons in the decay chain), 
but the significance to this shape fitting will diminish 
greatly without changes to the manner in which the gau- 
gino signal is extracted. This will necessarily degrade the 
accuracy with which mass differences between low-lying 
superpartners and the LSP can be determined. 



VI. CONCLUSIONS 

In this work we have provided a case study on how 
"New Physics" might appear in the early stages of 
the LHC experiment at CERN. If this new physics is 
the minimal supersymmetric standard model (or even 
more specifically, one of the specially-defined benchmark 
points of this model that have been designed by theo- 
rists), then it is very likely that the presence of the new 
physics and measurement of many key properties will 
be straightforward. It is not at all unlikely, however, 
that a supersymmetric world may be far from minimal. 
Evidence in support of this comes both from top-down 
constructions of GUT or string models as well as con- 
siderations of a bottom-up nature which focus on natu- 
ralness oftheelectroweak symmetry-breaking sector in 

SUSY [iUGME^. 

We have presented here just such an alternative world 
that is very much a plausible outcome of semi-realistic 
string constructions. The common element of these 
models is the presence of S'[/(3)-charged squarks and 
quarks which are isosinglets under weak SU(2) and 
carry equal and opposite hypercharge. Such states may 
present themselves as new 'partners' of Standard Model 
down-type quarks and squarks, as leptoquarks or di- 
quarks, or may have more exotic properties which for- 
bid renormalizable-level couplings altogether. We have 
used the language of Eq to characterize these states for 
convenience, but Eq GUT models are by no means the 
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only way in which such states can arise. The potential 
phenomenology of these states is both rich and tightly 
constrained, providing an interesting scenario for untan- 
gling high-scale physics from the low-energy observables 
of the LHC. 

If such states exist at LHC-accessible energies their 
presence will most surely be detected. This was the 
conclusion of previous work that considered such states 
in isolation, with only Standard Model fields serving as 
'background.' In our work we have provided the first 
analysis of such fields in the simultaneous presence of 
'standard' superpartners, using the benchmark point of 
SPS la as our example of simple SUSY models. The re- 
sult of our study has been to point out important conse- 
quences for experimenters hoping to observe and measure 
the properties of new physics in both sectors. A strong 
supersymmetric signal will not be obscured by these ex- 
otic states, but recognition of the new states may not be 
immediate. In the interim, we can expect fits to stan- 
dard SUSY-breaking paradigms such as minimal super- 
gravity, to fit well initially. As the LHC experiments inte- 
grate luminosity, we expect statistically significant devia- 
tions from the predictions of the pure MSSM to arise. In 
the leptoquark exotic case considered here that deviation 
will first appear in events with both leptons and high-p^ 
jets, but with little missing transverse energy. We an- 
ticipate that astute experimentalists will consider a lep- 
toquark interpretation of these 'anomalous' events, but 
most SUSY analysis techniques (particularly those which 
seek to isolate certain cascade decay chains and therefore 
measure certain mass differences) will need significant re- 
finement. We hope that the work presented here will spur 
others to provide a more sophisticated treatment of the 
quite-plausible 'what-if ' scenario we provide, and to pro- 
vide similar checks on the robustness of long-established 
analysis techniques. 
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APPENDIX A: ALTERNATIVE QUASI-STABLE 
MODEL 

As mentioned in Section IIIB1 stability on the time- 
scale of the detector at a hadron collider can be ac- 
complished by forbidding dimension-four operators that 
couple the exotic SU(3) representations to the Standard 
Model. In Sections III Bl and IIII CI this was illustrated 
by an ^-motivated model in which the dimension-four 
operators were forbidden by a combination of B, L, and 
U(1)n symmetries. Here, we consider a toy model, which 





Ov 
Y 


n' 


r lclQ 




n' 


3 x Q 


1 /fi 


1 


H" 


1 12 


_2 


3 x u 


-2/3 


1 


Hrl 


-1/2 


_2 


3 x d c 


1/3 


l 


2xD t 


-1/3 


-3 


3 x L 


-1/2 


l 


2 x Dl 


1/3 


-3 


3 x e+ 


1 


l 


2 x L s 


-1/2 


-2 


3xu c 





l 


2xL 3 


1/2 


-2 


S 





4 


T 





6 



TABLE VIII: Particle Spectrum for the alternative 
Quasi-Stable Exotic Model. The fields of the quasi-stable 
toy model are listed with hypercharge U(1)y and exotic U(l)' 
charge assignments. The new fields L s and Dt come in two 
vector-like pairs. Dt transforms as a charge —1/3 color triplet, 
while L s and L c 3 can be interpreted as additional Higgs pairs 
or as exotic lepton doublets. The new fields S and T are 
Standard Model singlets. 



was first constructed in |14[, in which this is achieved 
through a new U(l)' gauge symmetry alone. We will im- 
pose U{1)' charge assignments such that our exotics are 
stable at the renormalizable level but are allowed to decay 
through higher-dimensional operators. To achieve this 
we introduce an anomaly-free U{1)' with chiral matter 
representations which include two 5 + 5 representations 
of SU (5) and various fields which will be singlets under 
the Standard Model (SM) gauge group. Since B and L 
are not imposed, we allow for the additional constraints 
imposed by the non-observation of proton decay. 

The field content of the model is given in Table IVIIII 
The J7(l)' charge assignments are obtained by solving the 
anomaly cancellation equations. For the U(l)' charges 
listed below all the mixed anomalies of G SM x U(l)' and 
Y x C/(l)' 2 vanish. Since we only include additional rep- 
resentations 5 + 5, we will not modify the SM gauge ano- 
maly cancelations. However, additional SM singlets are 
still needed to cancel the trace anomaly and the t/(l)' 3 
anomaly. There is some freedom in choosing their U(l)' 
charge assignments [l4|. Note that the model maintains 
the unification of the SM gauge couplings at one loop. In 
fact, complete unification can be achieved for the appro- 
priate rescaling of the U(l)' gauge coupling g'. However, 
there is no simple explanation of the U(l)' charge quanti- 
zation, analogous to the U(1)y charge quantization in or- 
dinary grand unification theories, since the U(l)' charge 
assignments are just numerical solutions to the anomaly 
cancelation equations. 

Aside from the usual Yukawa interactions of the 
MSSM, the superpotential terms allowed by the choice 
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of Q' charges in Table IVIIII are given by 



dim4 



dim5 



SH u H d , TD t D%, SL S L C S , 
H d Qd c , H u Qu c , H d Le+, H u Lv c , 
+(H U -» L c s ) + (H d - L s ), (A.l) 
D^QH d S, D c t QL s S, D c t QQu c , 
D t d c v c v c , D t u c e + v c , D c t QLv c , 
D c t u c u c e + , D t QQv c , D c t d c u c v c . (A.2) 



the proton decay constraint is satisfied provided that 
M > 10 15 - 10 16 GeV. Alternatively, if D c t QH d S is re- 
placed by D c t QQu c , then M > 10 12 - 10 13 GeV is suffi- 
cient. 

All of the dimension-five operators in (|A.2p allow the 
D t to decay. However, up to counting factors the results 
are expected to be similar to those for the i?6-motivated 
model discussed in Sections [II Bl and I lII CI 



The exotic particles D t and D c t are stable at renormaliz- 
able level: there is only one operator TDtD^ involving D t 
that is allowed by the symmetry, and one can therefore 
assign a conserved D t number. This operator can provide 
a mass term for the exotic quarks should a [/(l)'-breaking 
vacuum expectation value (T) develop. The L s and L c s 
doublets can decay by the dimension 4 operators. 

There are a number of D t -number violating interac- 
tions allowed at dimension 5, however. These operators 
also violate baryon and lepton numbers, leading to pro- 
ton decay. The operators can be divided into two sets. 
The first two rows of (|A.2|) would allow B conservation 
for the assignment B = |(— -|) for D t (D$), while the last 
row of (|A.2|) would require B = -§(+§)• Clearly, if op- 
erators from both sets are present B will be violated and 
the proton will decay. However, since the _D t -number is 
conserved at dimension-4 level (and the D t is assumed 
to be heavier than the proton) , the proton decay rate T p 
is proportional to M~ 4 , where 1/M* is the scale multi- 
plying the dimension-5 operators. Depending on which 
operators are included, the dominant decay processes are 
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+v c (A B -l =0), or the 
corresponding neutron decays n — > ir°i/(i/ c ). 19 The ex- 
perimental limits for these process es are r > 10 31 yr (p — > 
n+v) and r > 10 32 yr (n — > n v) jl37l |. If the operator 
DtQH d S is present, with the scalar components of H d 
and S acquiring VEVs, as well as D t QQv c or D^d c u c v c , 

19 Decays into e+ are highly suppressed because the operator 
D^u c u c e + must be antisymmetric in the u c u c flavor indices, 



APPENDIX B: EXPLICIT CROSS-SECTION AND 
PARTIAL WIDTH EXPRESSIONS 

In keeping with the conventions of PYTHIA wc quote 
parton-level cross sections for 2^2 processes as differ- 
ential cross-sections of the form da/dt. In what follows 
the kinematic factor ^34 is given by 
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where the last relation holds when 7713 = = Mb- Wc 
will also have occasion to use the quantity 
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In PYTHIA the parton-level Mandelstam variables are 
given by 
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where 9 is the polar angle of parton 3 in the cm. frame 
of the hard scattering. 

leading to decays including a virtual c or c c quark. 



Production Cross Sections 



Q —* D 1/2 D 1/2 



There are two parton-level sub-processes of relevance: the q-q annihilation (two quarks of the same flavor) 
through gluon exchange in the direct channel and a diagram with t-channel exchange of a squark or slepton 
through the Yukawa interaction. The differential cross-section for the leptoquark case is given by 
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where a y = A 2 /4-7T. For the diquark case we have 
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Note the following: (1) In doing the PYTHIA analysis we always take A 6 = A 7 , A 8 = and A 9 = A 10 . So 
expressions involving a y should be understood as involving the common value of this Yukawa interaction. (2) 
For the Snowmass point we use as our superpartner spectrum, the up and down type left-handed squarks have 
the same masses to within a percent and up and down right-handed squarks have the same masses to within less 
than a percent. Thus, in evaluating (|B.5[) PYTHIA takes just the (common) gx soft mass and q~R soft masses. 

g + g -> D 1/2 D 1 /2 

There are two sub-diagrams for this process: the diagram with the initial gluons leading to a gluon intermediate 
state and then to two fermions, and the t-channel/u-channel exchange of the heavy fermion -Di/2- The differential 
cross-section is given by 
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plus the equivalent expression with i u. 
q + q —> DqD q 

For this process there are two types of diagrams: q-q annihilation (two quarks of the same flavor) through gluon 
exchange in the direct channel and t-channel exchange of a Standard Model fermion. As we allow our exotics 
to interact solely with first generation Standard Model fermions, the appropriate differential cross-section is 
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when q,q are a flavor other than the first generation and 

da _ da ira s a y J (M| 4 - t)(ii - i) + s(M| 4 + t) 
~dt ~ ~df + 18s 2 \ §t 
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for first-generation q,q. 
g + g^D Q + Do 

The QCD diagrams for the diquarks and leptoquarks are the same. There are several sub-diagrams: the diagram 
with the initial gluons leading to a gluon intermediate state and then to two scalars, the t-channel/u-channel 
exchange of the heavy scalar Do, and the "sea-gull" four-point diagram. The appropriate differential cross- 
section is given by 
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q + g -> D + f 

Single scalar leptoquarks and diquarks can be produced in association with a lepton and quark, respectively. 
The diagram involves gq initial states, either a quark in the direct channel or an exotic squark in the indirect 
(u) channel, and then a scalar and fermion in the final state. One of the two vertices is then given by QCD 
while the other is from the appropriate Yukawa coupling. The expression for the leptoquark is given by 
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The equivalent expression for the diquark gets an extra factor of two from the different color contraction 
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q + q 



This case is slightly different. What's coded into PYTHIA is (s/n)da, and this is given in turn by 
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2. Two-Body Decays 

For two-body decays of the form a — > & + c let us define the final momentum in the rest frame of a as 

[m 2 a - (m& + m c ) 2 ] 1/2 K - (™6 - ™c) 2 ] 1/2 
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Furthermore, we define the exotic scalar mixing matrix by the relations 

D = Dlcos9 D -Dlsm6 D ; D" = sintfzj + D% cos0 D . 
Then the partial widths for the exotic decay channels considered are given by the following: 20 
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where the | is from color. 
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where Q^j = —4, 0:2 = f 2 -, and iV,i(7Vj2) is the mixing element relating %° to the bino (wino). We ignore 
contributions to decays from neutralinos associated with the Standard Model singlet S (singlinos) or additional 
Z 1 (Z'-gauginos). 
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For the process — » //', neglecting fermion A_dq, we have 

masses and taking X 6 = X 7 = \lq and A 9 = A 10 = 
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We always give the spin-averaged decay rate for the exotic fcrmi- 
ons. For QCD production, this ignores spin correlations between 
the produced D and D c . It also neglects possible polarizations 
for leptoquark or diquark production processes. 



T s = ^M Do (B.17) 
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<k = 2(1 + 3 cos 2 0£>)|Adq| 



(B.18) 



• When kinematically allowed, the rate for the pro- 
cess Di/2 — » //' is given in the same approximation 
above by 



1 - 



m), - 



M 



D 



1/2 J 



(B.19) 



with di = |Alq| 2 for all Dy 2 decays and d% — 

8|Adq| 2 , 2\Xdq\ 2 for D±% decays to left-handed 
and right-handed squarks, respectively. 



X* L = +^N 12 + i taxi9 w N n 



X UR = +- tan 9 w Nu 



X„ 



x* = 



1 



tan 0w Nu 



1 



N u - -t&.n6 w N n , 



(B.22) 
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3. Three-Body Decays 

Finally, we give the rates for three-body decays. Let us 
keep only the first generation of Standard Model fermi- 
ons and ignore Jl/r mixing, except for the exotic scalars 
whose mixing is defined in (|B.12[) . We continue to ignore 
contributions to decays from singlinos or Z'-gauginos. 
We will also ignore Higgsino contributions. We define 
the quantities 



**D L 



x r . 



1 



(B.20) 



1 1 

X dL = -jj^ia + e tan6W7Vii 



where N\\ and N\2 are the entries of the (4x4) neutralino 
mass matrix. In practice we will take Nu — 1 and N12 = 
0. 

Let us begin by considering the process — > Xid c ' uC 
via the operator with coefficient A 10 in (|2.4p . Diagrams 
leading to this final state can involve the virtual Dq ^ 
state or virtual right-handed squarks. The amplitude- 
squared, summed over final colors and averaged over ini- 
tial spins, is given by 



\M\^4C\X w \'gi\A\\ 



(B.27) 



(B.21) where C — 2 is the color factor for this process and 



I-4I 2 = \Y dR \ 2 (p u c -p D )(p d c -p x ) + \Y UR \ 2 { Pd o -p D )(p u c -p x ) + (\Y Dl \ 2 + \Y Dr \ 2 ) { Pu c -p d c){p x -p D ) 
-R- e (Y dn Y UR )[(p u c -p D )(p d c -P x ) + (Pu° ■Px)(Pd"- -Pd) - (p u ° -Pd°)(Px "Pd)] 
-R- e (YD R Y dR )[(p u c ■p D ){p d c -P x ) + (Pu- -Pd-)(p x - Pd) - (Pu- -P x )(Pd- • Pd)} 
-R- e ( Y D R Yu R )[(p d <-- -Pd)(Pu- -P X ) + (Pw -Pd<=)(p x - Pd) - (p d - -P x )(Pw - Pd)} 
+ {2Re(Y*Y DR ) - Re(Y* L Y dR ) - Re(Y* L Y UR )} M D M x {p u . ■ Pd .) . 

The Yi are given by the expressions 



Y dn - 



X„ 



(p u ° - Pd) 2 ~ m \ 



Y — 

1 mr — 



X,, 



(Pd- - Pd) 2 -m 2 ~ 



(B.28) 
(B.29) 



Y, 



Dr. 



sin 2 8d 



cos 2 8d 



(Px -Pd) 2 - m D i (p x - Pd) 2 - m 2 2 



Y Dr = - Xd l cos On sin 0d 



(Px-Po) 2 -m D i (Px-Pd) 2 - 



v dI 



(B.30) 



All outgoing momenta are assumed to be physical mo- menta. In the rest frame of the initial state D^j^ we 
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therefore have the products 

(Pu° -Pd) = M D E U 

(Pd<= - Pd) = M D E d 



(B.31) 



(B.32) 



( Px ■ PD ) = M D E X = M% - M D {E d + E u ) (B.33) 



(iV ■Pd') = - -Ml + M D (E U + E d ) (B.34) 



Z Dr = +Xd r cos9dsul9d x 
1 



(Px~Pd) 2 - m Dl 

1 



(Px-Pd) 



(B.43) 



For leptoquark processes we distinguish between rates 
to u e and d v final states. In analogy to the Y and Z 
factors above, we define the quantities 



Qe L — 



{Pu-Pd) 2 - mj 



(B.44) 



ij>uo-p x ) = -Ml--Ml-M D E d 



1,1, 



' Px) = - 2 M x ~ MdEu ' 



(B.35) 



(B.36) 



The decay rate can be computed by integrating the par- 
tial width 



dr 



Mp 
256tt : 



■\M\ 2 dx d dx v 



(B.37) 



where X4 = 2Ei/Mu. The numerical integration can be 
performed for m u = m d = ( but retaining m x ^ 0) by 
following standard techniques [lldj | . 

For massless quarks, there is no interference between 
the diagrams involving the A 10 vertex and those from the 
A 9 vertex. The latter may be obtained from the above 
by the substitutions 



|A 10 | 2 - 

COS#D <- 

x lL - 

Thus, for example, we have 



1 2 



4|A , 
sin 9 e> 



-X 



in 



(B.38) 



Rdr 



{p e - Pd) 2 - m( 



X* 



(Pd-PD) 2 ~ mj 



X 



dh 



S UR — 



(Pv - Pd) 2 - m 2 -^ 



{Pu-Pd) 2 - mj 



{Pe-PD? -ml 



(B.45) 
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(B.48) 
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The appropriate three-body decay rates can then be 
found following the above description using the substi- 
tution rules summarized in Table HXl 



Yd r — ► Z Dl , Y Dl 



Y, 



Z d L , 



Y, 



• Z D„ 



(B.39) 



where 



Ml 



X dL 



(Pw - Pd) 2 - mi 



(B.40) 



X„ 



{Pd- - Pd) 2 -mi 



(B.41) 



-Dr 



-X 



D, 



COS Ve> 



sin Ve> 



(Px-Pd) 2 - m D i (Px-Pd) 



(B.42) 
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Tio 


r s 


IT 


J- 7 


r 6 


D -> u c d c x 


D -> u c d c x 


D — > ue.\ 


D — ► d^x -D — * «ex 


(7IA 10 I 2 


4C|A 9 | 2 


IA 7 I 2 


\\ 7 \ 2 


X 6 I 2 


Pd c 


Pd c 


T>*> 
l'<- 




Pe 




Pu c 


Pu 


Pd 


Pu 


Y , 
^d R 


Y , 


Y* 


-x* L 


x i R 


** R 
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-XL 


X UR 


Xd l 


— Xd b 


X Dl 


Xd l 


—Xd r 


Xd r 


—Xd l 


Xd r 


Xd r 


—Xb l 


cos 6d 


sin 0_d 


cos 8d 


cos 8d 


sin Od 


sin Od 


cos 8d 


sin 6d 


sin #d 


COS #£) 


m\ 

or 




m% L 




m| H 






ml L 


771 j 

dh 


ml R 




z Dr 


Yn L 


Yd l 


Z D R 




z Dl 


Yd r 


y Dr 


Zd l 




Zd L 


Qe L 


R VL 




v 

u R 




Qu L 


Rd L 


S UR 



TABLE IX: Substitution rules for obtaining three-body de- 
cay rates. For each process, replace the quantity in JB . 271 ) - (IB.36I 
with the appropriate variable from the table. 
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